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Abstract 


A  new  interpretation  of  the  anomalous  dispersion  process  which  is 
observed  at  low  frequencies  is  presented  which  is  based  on  a  realistic 
description  of  the  structural  ordering  and  fluctuation  in  carrier  dominated 
dielectrics.  It  is. shown  that  this  form  of  response  occurs  for  systems  of 
low  spatial  dimensionality  and  generates  a  sample-size  dependent  conductivity 
The  relationship  of  the  mechanism  to  that  of  power  law  noise  in  electrical 
systems  is  identified  and  the  structural  interpretation  explored.  Particular 
features  of  hydrogen  bonded  systems  are  described  in  which  the  dispersion  is 
likely  to  be  important  in  a  biological  context.  It  is  pointed  out  that  the 
process  is  directly  relevant  to  solid-state  batteries  such  as  those  based  on 
fast-ion  conductors  such  as  6 -alumina. 


'I 


y 


Vr 


Accession  For 


NTI5  CR1*I 
DTIC  t;.s 
U:tarutoiui'.’sd 
Just  IF  lcatlo: 


Br- 


Dlatrlbutlon/ 


l - , 


Diet 


ft 


Availability  Codes 


Avail  and/or 


Spec lal 


1.  Introduction 


In  the  previous  report  (Dissado  and  Hill,  1982) the  effect  of  sample 
size  in  the  macroscopic  range  was  investigated  in  ferroelectric  type  materials 
as  specific  examples  of  bound  charge,  dielectrically  active,  systems.  In  the 
present  study  the  complementary  case  of  the  effects  of  sample  size  in  almost 
free  carrier  systems  have  been  investigated.  It  is  shown  that  the  anomalous 
low  frequency  dispersion  which  is  invariably  observed  in  quasi-free  ion 
carrier  systems  such  as  the  Hollandites  and  8-alumina  (Jonscher,  1978),  does 
show  characteristic  size  effects  which  can  appear  as  a  sample  size  dependent 
quasi-D.C.  conductivity.  Although  these  effects  have  only  been  characterised 
for  ionic  carrieu  it  is  argued  that  similiar  effects  can  occur  in  low 
dimensionality  systems  for  electron  carriers  but  will  occur  on  much  shorter 
time  scales. 

The  dielectric  response  of  a  solid  can  be  characterised  by  the 
dependence  of  the  complex  susceptibility,  x(“)  *  x'Oo)  “  ix"(<*0,  the 
frequency,  u,  and  on  the  external  variables  such  as  temperature,  pressure, 
humidity,  etc.  In  general  a  peak  is  observed  in  the  frequency  dependence  of 
the  loss  component  x"(“>)  together  with  a  dielectric  increment  in  the  real 
component,  x’  (<*>)»  corresponding  to  a  relaxation  process  (McCrum  et.al.,  1967). 
Typically  such  loss  peaks  have  a  non-Debye  shape  (Jonscher,  1975,  1977;  Hill, 
1981a)  which  is  given  by 

X '  (id)  «  x"(<»i)  *  «  ^"”n^  ;  u>  >  u)p  (la) 

and  x'^w)  *  x(°)  “  a.x'(u)  «  w*™  ;  to  <  u>p  (lb) 

where  a  is  a  constant  and  x(0)  is  the  susceptibility  at  zero  frequency.  The 
constants  m  and  n  lie  between  zero  and  unity  and  u  is  the  relaxation  rate 
and  may  be  identified  with  the  frequency  of  maximuin  loss.  A  consequence 
of  the  high  frequency  behaviour  of  equation  (la)  is  that  the  ratio  x"(«)/x'  (<•>) 
is  independent  of  frequency  (Cole  and  Cole,  1941;  Jonscher,  1975)  and  is 
given  by 

X"(«)/X'(u)  -  cot(nir/2)  (2) 

However  Jonscher  (1978)  has  identified  a  second  class  of  behaviour  in  which 
the  frequency  dependence  of  equation  (la)  obtains  through  a  transition  from 
a  small  value  of  the  index  n  at  low  frequencies  to  a  different,  and  larger, 
value  beyond  a  characteristic  frequency  u  ,  as  shown  in  figure  1.  Relationship 
(2)  is  obeyed  in  both  the  high  and  low  frequency  limits  with  the  appropriate 
value  of  the  index  n.  Neither  frequency  range  has  yet  been  observed  to 
form  part  of  a  loss  peak,  furthermore  x(«)  has  been  shown  to  be  normalisable 
(Hill,  1978)  with  respect  to  variations  in  u  due  to  changes  in  temperature 
(Jonscher,  1978)  and  measuring  field  (Jonscher,  Meca  and  Millany,  1979)  in 
a  similiar  manner  to  a  loss  characteristic  (Ferry,  1970;  Hill,  1978).  The 
conditions  for  this  to  be  possible  should  u  simply  represent  the  point  of 
overlap  of  two  characteristics  are  very  restrictive,  namely  that  the 
susceptibility  increment  and  characteristic  frequency  for  each  process  should 
behave  in  an  identical  fashion  with  respect  to  variations  in  temperature  and 
field,  while  yet  being  different  in  magnitude. 


(2) 


Thus  this  type  of  behaviour  can  best  be  regarded  as  a  feature  in  its 
own  right,  and  since  it  is  the  lowest  frequency  process  observeable 
dielectrically  it  has  been  termed  the  anomalous  low  frequency  dispersion 
by  Jonscher  (1978).  For  convenience  we  define  the  value  of  (1-n)  in  the 
lower  frequency  range  to  be  p,  and  observe  from  the  literature  quoted  that 
p  lies  in  the  range  0.7  to  1.0. 

When  avai  lable  frequency  range  only  extends  down  to  the  dotted 
line  in  figure  1  the  phenomena  may  be  confused  with  the  onset  of  D.C. 
conduction,  for  which 


X  *  (oj)  =  constant 
x"(u))  «  Ojj  c/li! 


u  0 


(3a) 

(3b) 


with  oQC  the  frequency  independent  D.C.  conductivity.  This  behaviour  is 
indicated  by  the  light  dotted  curves  in  figure  1.  Whilst  the  presence  of 
a  D.C.  component  cannot  be  ruled  out  the  validity  of  relationship  (2)  in 


the  extremely  low  frequency  region  effectively  negates  this  as  an  explanation 
of  the  dispersion  process.  Similiarly  the  Maxwell-Wagner  effect  (Maxwell, 
1954;  Wagner,  1914)  predicts  a  limiting  behaviour  of  the  form  of  relationship 
(3)  but  with  a  higher  frequency  region  in  which 


x"(<*>) 


X'M 


with  the  ratio  of  the  imaginary  and  real  components  proportional  to  the 
frequency. 


It  has  been  proposed  (Poliak  and  Geballe,  1961)  that  the  high  frequency 
portion  of  the  anomalous  dispersion  can  be  described  in  terms  of  variable 
range  hopping  of  carriers  with  values  of  n  estimated  to  lie  between  0.5  and 
1.0  (Austin  and  Mott,  1969;  Hill  1977).  This  model  however  offers  no 
convincing  explanation  of  the  lower  frequency  dispersion  which  is  a  major 
feature  of  the  response.  A  modified  version  of  the  model  predicts  a 
temperature  dependence  of  n  (Elliot,  1977)  which  is  rarely,  if  ever,  observed 
(Frost  and  Jonscher,  1975;  Hill  1981a).  Clearly  such  an  explanation  is  partial 
at  best.  Jonscher  (1977)  has  suggested  that  both  regions  of  the  anomalous 
response  should  be  regarded  as  examples  of  the  power  law  response  common  to 
all  dielectrics,  for  which  an  explanation  was  suggested  by  Ngai  et.al.  (1979) 
in  terms  of  a  cooperative,  many-body ,  model.  It  was  suggested  that  the  two 
regions  of  behaviour  could  be  regarded  as  resulting  from  independent  processes, 
although  the  nature  of  the  difference  was  not  specified.  Furthermore  in  a 
few  cases  (Jonscher,  1978)  the  power  law  of  the  anomalous  low  frequency 
dispersion  undergoes  a  further  modification  at  the  lowest  frequencies  to  a 
characteristic  with  a  value  of  p  of  one  half.  This  latter  fact  cannot  be 
explained  on  the  basis  of  a  third  independent  process,  since  such  a  process 
would  not  be  observeable  in  the  presence  of  the  dispersion  with  p  close  to 
unity.  Such  behaviour,  where  it  occurs,  must  be  regarded  as  the  extremely 
low  limit  of  the  anomalous  dispersion,  and  the  power  law  exponent  of  0..5 
strongly  suggests  that  it  is  determined  by  independent  particle  diffusion. 

The  full  explanation  of  this  class  of  dielectric  response  therefore 
seems  likely  to  involve  a  composite  process  whose  component  parts  can  be 
identified  with  dynamical  possibilities  inherent  in  the  material  structure, 
in  a  similiar  manner  to  that  of  the  loss  peak  (Dissado  and  Hill,  1979;  Dissado, 
1982),  and  possessing  the  observed  limiting  particle  diffusion  behaviour.  We 
propose  here  that  these  criteria  can  be  met  by  a  model  based  on  the  concept  of 


a  cluster  of  quasi-mobile  charges  possessing  partial  structural  regularity, 
for  which  the  model  of  Ngai  et.al.  (1979)  is  a  particular  case  of  a  more 
general  description  which  leads  to  a  power  law  response  (Halperin  et.al., 
1972;  Binder  et.al.,  1975;  Patashinskii  and  Pokrovskii,  1979).  In  the 
following  sections  we  intend  to  identify  the  necessary  elements  for  this 
model  in  the  systems  that  produce  this  class  of  response  and  delineate  the 
maimer  in  which  they  contribute  to  the  combined  behaviour.  During  this 
exposition  the  connection  with  other  electrical  phenomena  such  as  electrical 
noise  (Gupta,  1977;  Hill  et.al.,  1981)  and  D.C.  conduction  will  become 
apparent. 


2.  Conditions  necessary  for  Low  Frequency  Dispersion 

The  class  of  response  identified  above  as  anomalous  low  frequency 
dispersion  has  been  observed  in  a  wide  range  of  materials  of  considerably 
different  chemical  nature.  These  include  fast  ion  conductors  such  as 
Hollandite  (Jonscher,  1978),  ionic  glasses  (Jonscher  and  Frost,  1976;  Doyle, 
1981),  wetted  proteins  (Jonscher  1978;  Eden  et.al.,  1980),  polymers 
(Bodiakan,  1982),  solid  nickel  suspension  (Ramdeen,  1981),  humid  sand 
(Shahidi  et.al.,  1975)  and  stearic  acid  films  (Jonscher  et.al.,  1979).  A 
comparison  between  the  response  of  humid  sand  and  sand/water  mixture  revealed 
that  the  dispersion  is  associated  with  the  water  present  as  an  adsorbate. 

Since  the  catalytic  activity  of  the  solid  nickel  suspension  depended  upon 
the  adsorption  of  molecules  it  is  likely  that  in  this  case  the  effect  also 
relies  on  the  presence  of  water  as  an  adsorbate.  Furthermore  the  observation 
of  the  low  frequency  response  in  hydrated  proteins  has  been  strongly  linked 
with  the  filling  of  secondary  hydration  sites  (Eden  et.al.,  1980).  In  the 
case  of  polymers  it  is  known  that  water  diffuses  into  materials  relatively 
easily  (Auckland  and  Cooper,  1974)  and  adsorption  may  occur  at  the  interface 
between  the  crystalline  and  amorphous  regions  such  as  spherulite  and  lamella 
surfaces  (Barham  and  Keller,  1977).  It  therefore  seems  probable  that  the 
observation  of  the  dispersion  in  all  hydrated  samples  is  a  consequence  of 
the  formation  of  adsorbed  water  layers  on  the  interior  surfaces  of  one  of 
the  components  of  a  heterogeneous  material  (Jonscher,  1978). 

Although  ionic  materials  such  as  Hollandite  possess  a  distinctly 
different  structural  nature  from  that  of  the  hydrated  systems  there  are  some 
significant  similiarities.  In  the  Hollandites  the  crystal  structure  forms 
channels  through  which  a  migrating  ion  may  pass  to  give  D.C.  conduction. 

These  channels  £re  one-dimensional  in  the  Hollandites  and  are  partially 
occupied  with  K  ions  (Reau  et.al.,  1977;  Jonscher  et.al.,  1979).  The 
surrounding  crystal  structure  creates  a  sinusoidal  potential  surface  along 
the  channels  with  minima  at  regular  separations  of  2.9  A  and  a  potential 
barrier  of  approximately  0.38  eV  (Beyeler,  1976).  In  each  crystal  unit  cell 
there  are  four  channels  in  parallel,  figure  (2a),  with  a  correlation  length 
along  the  channels  of  approximately  35  A,  composed  of  channel  unit  cell  groups 
of  3  to  6  minima  each  containing  a  vacancy,  figure  (2b),  (Beyeler,  1976). 
Because  of  the  strong  ion-ion  interaction  the  regularity  of  the  spacing  is 
distorted  with  an  estimated  displacement  of  about  22%  neighbouring  the  vacancy. 
These  supercell  groups  are  combined  together  so  as  to  satisfy  the  material 
stochiometry,  giving  a  gross  structure  which  can  only  be  described  in  terms 
of  a  probability  distribution  of  supercell  types.  The  similiarity  with  the 
adsorbate  system  is  now  obvious,  since  the  characteristic  features  of  these 
are  the  presence  of  a  sinusoidal  surface  potential,  produced  by  the  fixed 
structure  of  the  surface,  and  a  partial  occupancy  of  the  minima  by  the  adsorbed 
species.  Furthermore  when  the  occupancy  is  of  the  order  of  0.5  it  would  be 
reasoneable  to  expect  intermolecular  interactions  such  as  hydrogen-bonding 
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to  distort  the  local  structure  of  the  surface  potential  which,  however,  is 
sufficiently  strong  to  impose  a  near  regular  molecular  ordering. 

When  we  examine  the  other  materials  that  exhibit  the  anomalous  type  of 
behaviour  it  is  found  that  they  fit  into  the  same  pattern  of  system.  The 
ionic  glasses  are  heterogeneous  materials  with  crystallites  providing  an 
interior  surface  to  the  amorphous  regions,  whereas  the  stearic  acid  fims  are 
surface  adsorbed  via  the  organic  acid  (  -CO_H  )  groups  which  are  likely  to 
have  a  partial  occupancy,  particularly  if  charge  injection  from  the  electrode 
ionise  some  of  the  groups. 


3.  The  Mechanism  of  Low  Frequency  Dispersion 

In  the  preceeding  section  we  have  identified  the  low  frequency  dispersion 
response  with  systems  in  which  a  mechanically  rigid  lattice  forms  a  regular 
array  of  binding  sites  for  ions  or  molecules,  most  of  which  are  occupied.  This 
large  occupancy  factor  allows  the  bound  species  to  interact  with  one  another 
so  as  to  modify  the  regularity  of  the  inter-site  spacing  from  that  of  the 
substrate  potential.  As  a  result  the  translational  reproduceability  will  be 
destroyed  over  a  distance  £  which  defines  a  positional  correlation  length 
for  the  cluster  of  ions  or  molecules  thus  formed.  For  separations  less  than 
£  patterns  of  local  ordering  will  exist,  such  as  the  supercell  groups  in  the 
HSllandites.  In  the  water  adsorbate  systems  it  is  likely  that  the  local 
ordering  is  that  of  hydrogen  bonding  molecule  groups  whose  network  regularity 
is  over-ridden  at  £c  by  the  greater  spacing  of  the  adsorbate  sites. 

If  substrate  systems  such  as  these  were  sparsely  occupied  by  a  low 
concentration  of  ions  charge  transport  could  occur  by  particle  hopping.  Since 
each  ion  would  have  unrestricted  access  to  an  unoccupied  site  on  a  path 
connecting  the  electrodes  a  D.C.  conductivity  would  be  expected  from  them. 

When  the  site  occupancy  is  high,  however,  access  to  sites  on  an  electrode 
connecting  pathway  is  substantially  reduced  when  the  dimensionality  of  the 
system  is  low.  As  a  result  charge  transport  in  the  highly  occupied,  low 
dimensionality,  systems  must  be  considered  in  an  entirely  different  manner  to 
that  of  the  unrestricted  hopping  systems.  The  approach  adopted  here  is  a 
modification  of  that  previously  used  to  describe  relaxation  in  condensed 
matter  (Dissado  and  Hill,  1979,  1980,  1981;  Dissado,  1982). 

In  this  picture  materials  exhibiting  relaxation  phenomena  are  envisaged 
as  being  composed  of  clusters  in  which  the  perfect  lattice  structure  is 
distorted.  A  macroscopic  sample  of  the  material  will  be  formed  with  a  steady- 
state  distribution  of  such  clusters  varying  continuously  in  form  from 
complete  distortion  to  a  perfect  lattice  structure.  This  distribution  should 
be  viewed  in  the  dynamic  sense  in  which  a  particular  cluster  if  labelled  would 
be  seen  to  adopt  all  the  possible  structures  with  a  fixed  probability  density 
over  a  suitable  time  scale.  This  description  should  not  be  confused  with 
that  of  a  static  distribution  in  which  the  structure  would  vary  from  site  to 
site  but  remain  constant  in  time.  The  mechanism  whereby  individual  clusters 
fluctuate  their  structure  is  via  inter-cluster  (I-C)  exchanges,  which  include 
the  transport  of  distortion  or  strain  by  way  of  a  virtual  phonon  process 
(Joffrin  and  Levelut,  1975)  in  addition  to  any  particle  exchanges  that  may  be 
possible. 

The  response  of  such  systems  to  an  externally  applied  spatially  uniform 
field  is  determined  by  the  relaxation  of  the  spatially  uniform  polarisation 
fluctuations  (Kubo,  1957).  This  relaxation  is  a  composite  process  in  which 
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the  evolution  of  the  cluster  structure  from  the  spatially  uniform  displacement 
is  initially  dominant,  to  be  followed  at  longer  times  by  the  evolution  of  the 
steady  state  cluster  distribution  (Dissado,  1982). 

We  have  previously  considered  the  situation  in  which  the  I-C  exchanges 
do  not  themselves  constitute  a  linear  response  of  the  applied  field.  As  a 
result  they  act  not  as  a  noise  (polarisation  fluctuation)  on  the  material 
polarisation  but  as  a  noise  on  the  relaxation  current  of  the  polarisation 
(Hill  et.al.,  1981),  retarding  the  approach  to  the  dynamic  equilibrium  and 
thereby  producing  the  observed  loss  peak  response  (Jonscher  et.al.,  1980; 
Dissado,  1982).  However  when  the  I-C  exchanges  carry  a  charge  such  as  in  the 
systems  considered  here  they  constitute  an  element  of  electrical  current 
which  couples  linearly  with  the  applied  field.  It  is  this  difference  which 
is  the  essential  characteristic  of  the  anomalous  low  frequency  dispersion 
response  which  is  observed. 

3.1  Inter-cluster  exchanges  and  the  effective  transport  of  charge 

Because  of  the  cluster  nature  of  these  systems  the  effective  hopping 
of  an  ion  between  available  sites  will  have  a  different  effect  depending  on 
whether  the  distance  covered,  L,  is  less  than  or  greater  than  the  cluster 
correlation  length  £  .  When  L  is  less  than  £  the  ion  will  still  be  partially 
correlated  with  its  original  site  containing  £he  counter  charge  and  these 
processes  must  be  regarded  as  polarising  the  cluster.  Only  when  L  is  greater 
than  £  do  the  ion  and  counter  charge  reside  in  different  clusters  and  move 
independently  of  each  other.  This  process,  which  corresponds  to  cluster 
ionisation,  generates  in  addition  a  new  internal  structure  for  the  clusters 
evolved.  In  the  Hollandites,  for  example,  the  supercell  composition  of  the 
clusters  will  be  different  immediately  following  the  interchange  although  an 
internal  readjustment  of  the  clusters  to  the  presence  of  a  hole  or  an  extra 
ion  will  result  from  the  ion-ion  interactions.  In  its  absence  the  inter¬ 
cluster  exchange  process  constitutes  an  electrical  current  element  which 
simultaneously  creates  a  range  of  different  cluster  structures.  Therefore 
a  description  of  the  evolution  of  the  probability  density  of  cluster  structures 
will  also  give  the  distribution  of  the  current  elements  produced  by  the  charge 
carrying  I-C  exchange  processes.  This  is  achieved  by  determining  the 
probability  of  realisation  of  a  prepared  initial  state  of  the  system  at  a 
subsequent  time  and  convoluting  it  with  the  probability  of  realisation  of  a 
new  state  of  the  system  (Hill  et.al.,  1981;  Dissado,  1982). 

The  convolution  is  necessary  because  in  a  closed  or  conserved  parameter 
system  decay  of  the  information  placed  in  the  initial  description  must 
inevitably  mean  the  growth  of  that  information  in  a  redistributed  form  within 
the  system,  no  mechanism  being  available  to  allow  loss  to  an  external  sink. 

In  this  case  the  initial  prepared  state  consists  of  a  translationally  equivalent 
set  of  identical  neutral  clusters  whose  internal  structure  satisfies 
stochiometric  requirements.  The  time  development  of  the  prepared  state  is 
determined  by  the  inter-cluster  exchange  processes  acting  as  a  perturbation  on 
the  system.  In  the  interaction  representation  the  time  development  is  given  by 

<Kt)  ■  Real  C  expE  “  exp{  F(t)}  )  (4) 

where  EE  is  the  average  energy  required  to  create  the  prepared  state  structure 
from  the  steady  state  equilibrium,  and  F(t)  describes  the  effect  of  inter¬ 
cluster  exchange  perturbations  V(t)  in  terms  of  their  cumulants  (Kubo,  1962; 
Chester,  1963) 
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rn  -iEt 

N(E){  1  -  e  lLt  >dE  (7) 
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where  E  *  E^  -  E^^  and  N(E)  is  the  density  of  energy  changes  Ee  produced  by 
the  perturbation.  The  form  of  N(E)  is  particularly  important  in  determining 
F(t)  and  can  be  deduced  in  the  present  case  by  considering  the  possible  ways 
in  which  a  given  amount  of  energy  can  be  utilised  in  producing  inter-cluster 
exchanges . 

The  transport  of  charge  between  clusters,  defined  as  an  I-C  exchange, 
possesses  a  number  of  possibilities  in  addition  to  the  direct  hopping  of  an 
ion  between  the  respective  clusters.  There  will  be  a  continuum  of  these 
ancillary  processes  in  which  correlated  I-C  exchanges  over  shorter  distances 
combine  to  give  an  effective  charge  transport  over  a  long  range.  The  number 
of  such  I-C  exchanges  required  to  produce  an  equivalent  transport  of  charge 
will  increase  as  their  transfer  range  decreases,  figure  3(a).  Since  the 
energy  required  to  separate  a  charge  from  its  counter  charge  increases  with 
separation  the  number  density  N(E)  of  inter-cluster  exchanges  at  a  particular 
energy  E  will  increase  as  E  decreases  and  the  leading  term  in  N(E)  will  have 
the  form 

N(E)  -  p/E  (8) 

which  leads  to  an  infra-red  divergent  behaviour  (Hopfield,  1970)  for  F(t) 

F(t)  -  ipnc  -  p{Ex(int)  +  ln(int)  +  y>  (9) 

The  average  energy  of  excitations  produced  by  the  perturbation  is  thus 


rn 

EE. 
•  n 


N(E)  dE 


and  since  Efj  is  the  maximum  energy  change  produced  by  the  perturbation,  that 
is  the  cluster  ionisation  energy,  the  index  p  must  be  unity  or  less. 

Each  of  the  I-C  exchanges  allowed  for  in  F(t)  alters  the  internal 
structure  of  the  clusters  involved  and  thus  the  factor  exp{F(t)}  in  equation 
(4)  describes  the  development  of  the  distribution  of  cluster  structures 
as  a  self-energy  clothing  the  prepared  state.  Because  the  prepared  state 
requires  the  input  of  the  average  energy  Ee  in  order  to  create  it  from  the 
distribution  it  is  this  energy  that  is  available  through  the  perturbations 
V(t)  for  the  I-C  exchanges  creating  the  new  structure  and  transferring  the 
charges,  hence 

EE  *  K  <E>  “  pEn  (11) 

The  energy  stored  in  the  ordered  cluster  will  be  present  as  delocalised 
vibrations  whose  quanta  are  in  excess  of  the  localised  vibrations  of  the 
disordered  cluster.  An  alternative  view  of  the  evolution  can  be  given  in 
terms  of  the  required  changes  in  the  vibrational  quanta  (Dissado,  1982). 


The  function  tKt)  should  be  regarded  in  the  quantum  mechanical  sense 
as  the  probability  of  realisation  of  the  prepared  state  at  time  t,  and  is 
presented  in  its  entirety  in  figure  A.  This  function  approaches  zero  time 
as  a  gaussian  with  time  reversal  symmetry  but  when  nt  *  1  partial  recurrence 
cycles  set  in  as  energy  is  exchanged  between  the  prepared  state  and  the  X-C 
exchange  fluctuations.  At  times  in  excess  of  nt  being  unity  a  power  law 
decay  dominates  as  the  large  number  of  different  I-C  exchange  processes 
involved  in  the  self-energy  overdamp  the  prepared  state  irreversibly  and 
clothe  it, 

i/>(t)  ^  exp{  -  (nt/2)2  >  ;  nt  <  1  (12a) 

*■  .  nt  >  1 
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(12b) 


It  should  be  noted  that  the  long  time  decay  of  equation  (12b)  has  been 
obtained  without  the  aid  of  an  arbitrarily  introduced  decay  factor  in  the 
correlations  of  V(t)  contributing  to  F(t)  as  is  usual  (Sesibois  and  De-Leener, 
1966)  in  problems  concerning  conservative  systems.  The  long  time  development 
of  iHt)  in  equation  (12b)  shows  that  the  probability  of  finding  the  prepared 
state  at  a  time  t  decays  as  the  alterations  in  cluster  structure  brought 
about  by_I-C  exchanges  grow.  The  probability  of  realising  a  clothed  state 
at  time  t  after  it  has  started  to  form  must  therefore  take  the  form 

_  _  P  _  Q3) 

♦  (t)  'v*  (nt)  ;  nt  >  1 

The  probability  density  of  cluster  structures  observed  at  time  t  can  be 
obtained  from  the  contributions  of  each  I-C  exchange  event  which  occur  with 
equal  facility  in  each  identical  time  interval.  The  contribution  from  an 
I-C  exchange  occurring  at  time  t.  is  given  by  the  probability  of  realising 
the  new  state  at  a  time  (t  -  t-)Hater  weighted  by  the  probability  of  finding 
the  prepared  state  at  time  t^.  This  contribution  is  thus 


P(t,t^)  —  .\«-  fcj# 

and  the  average  structure  is  therefore 

£p(t,t..) /Number  of  elements  -*■  t  *  [  t7^  (t  -  t.)*1*  dt..  (15) 

tx  Jo 

as  the  number  of  elements  is  proportional  to  t. 

The  establishment  of  the  steady  state  distribution  is  presented 
schematically  in  figure  5,  with  a  selected  number  of  contributions  shown. 

The  probability  density  of  cluster  structures  can  be  derived  from  equation 
(15)  (Hill  et.al.,  1981)  and  is  shown  in  figure  6  for  a  range  of  values  of 
the  index  p.  Although  P(t,t.)  only  allows  for  single  I-C  exchange  events 
the  probability  distribution  obtained  has  been  shown  to  be  persistent  when 
multiple  events  are  included  (Hill  et.al.,  1981)  and  thus  represents  a  steady 
state  description  of  the  system. 

The  fluctuations  producing  the  steady  state  distribution  are  driven  by 
exchanges  in  stored  energy  between  the  competing  features  determining  the 
short  range  and  long  range  structure  respectively.  In  the  case  of  the 
Hollandites  these  are  the  intra-cluster  ion-ion  repulsions  and  the  nett 
attraction  between  charged  clusters.  The  resulting  probability  density 
therefore  detetrmines  the  configuration  entropy  of  the  inter-cluster  super¬ 
structure  and  the  inter-cluster  exchanj  '8  produce  fluctuations  in  the  inter¬ 
cluster  configuration  entrop  •  Bees’  .  a  configurational  entropy  contribution 
to  the  probability  of  realisxnn  s,  ;em  is  temperature  independent  for  a  given 


’.(t  -  t,) 
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state  of  the  system  the  index  p  resulting  from  its  fluctuations  will  also  be 
temperature  independent. 

It  is  the  special  feature  of  the  systems  examined  here  that  the  I-C 
exchange  processes  transport  an  effective  electrical  charge  and  therefore 
constitute  an  elemental  electric  current.  Each  contribution  P(t,t^) 
therefore  represents  a  current  element  in  the  complete  distribution  as 
observed  at  time  t.  This  observed  current  can  therefore  be  obtained  by 
averaging  over  the  normalised  cluster  distribution 

<J>  -  Cr(l+p)  r(l-p)t3_1  f*  j  t7<t  -  t1)+p  dt.  -  j  (16) 

Jq  1  1  1 

where  j  is  the  current  element  which  is  given  by  the  rate  of  transport  of  the 

effective  cluster  charge  over  the  correlation  distance  £  .  This  expression 
can  be  best  understood  if  each  contribution  P(t,t^)  is  taken  to  represent  a 
number  of  sequential  nearest  neighbour  cluster  transports  of  effective  charge 

in  the  time  range  (t  -  tj) ,  the  nett  transport  must  then  be  divided  by  the 

number  of  primitive  transports  in  time  t,  to  give  the  average  current. 

3.2  The  internal  relaxation  of  clusters 

The  average  current  <J>  given  by  equation  (16)  has  the  form  of  a  D.C. 
current  with  fluctuations  which  give  rise  to  power  law  noise  (Hill  et.al.,  1981) 
which  originates  in  the  reorganisation  of  long  range  inter-cluster  structure. 

In  deriving  this  expression  we  have,  however,  neglected  the  possibility  of 
recombination  of  the  effective  charge  and  counter  charge.  Since  recombination 
can  only  take  place  for  separations  less  than  the  cluster  correlation  length 
£  the  dynamics  of  this  process  are  identical  to  the  relaxation  of  a  polarised 
cluster. 

The  time  development  of  the  relaxation  function  of  a  polarised  cluster, 
<f>(t)  has  already  been  obtained  (Oissado  and  Hill,  1979;  Dissado,  1982)  as 

<Kt)  'v  (Ct)-*  exp{  -  u>ct  }  (17) 

which  can  be  thought  of  as  the  possibility  of  realising  the  initial 
polarised  state  at  a  time  t  after  relaxation  has  been  initiated.  Here  u  is 
the  relaxation  rate  of  the  cluster  dipole  rather  than  a  site  dipole  and  §he 
(Ct)^1  factor  results  from  the  evolution  of  the  initial  state  into  the 
distorted  array  that  comprises  the  cluster  structure. 

In  the  system  considered  here  each  ion  within  a  neutral  cluster  can  be 
considered  to  form  a  site  dipole  with  its  lattice  counter  charge.  The 
polarisation  fluctuation  coupled  to  the  spatially  uniform  applied  field  consists 
of  completely  in-phase  displacements  of  the  site  dipoles  along  the  direction 
of  the  field.  This  is  equivalent  to  the  displacement  of  an  effective 
cluster  charge  over  a  distance  £  ,  figure  3(b),  just  as  the  I-C  exchanges  are 
equivalent  to  an  effective  displacement  of  this  cluster  charge  over  macroscopic 
distances. 

The  spatially  uniform  fluctuation  can  be  considered  to  be  a  prepared 
state  of  the  cluster  which  will  decay  into  the  distorted  structures  available. 
This  behaviour  can  be  described  as  in  the  previous  section  by  considering  the 
initial  state  to  be  created^  from  the  average  steady  state  cluster  structure 
with  an  amount  of  energy  RE  stored  in  the  structural  ordering.  This  energy  is 
available  through  perturbations  V(t)  to  reorganise  the  cluster  structure  into 
the  available  distorted  forms.  As  previously  there  will  be  an  increasing 
number  of  ion  displacement  processes  as  the  energy  requirement  for  the  individual 
process  becomes  less.  This  will  lead  to  an  infra-red  divergence  as  previously 
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discussed  and  give  the  time  power  law  decay,  (St)-11,  of  the  initial  state 
found  in  equation  (17).  Here  the  average  energy  of  fluctuation  formation 
is 

HE  *  nh£  with  0  (.n  i  1  (18) 

and  5  is  the  maximum  energy  change  required  by  a  single  cluster  distorting 
process.  In  the  present  case  £  is  the  energy  change  introduced  when  an  ion 
is  moved  into  a  neighbouring  vacancy  with  all  the  other  cluster  ions  undisplaced 
and  it  may  involve  details  of  the  local  potential  surface  as  well  as  ion-ion 
repulsions.  The  processes  involved  in  the  t  decay  are  those  establishing 
the  steady  state  configurational  entropy  density  of  neutral  cluster  structures 
and  thus  we  may  expect  n  to  be  also  temperature  independent  unless  a 
structural  ordering  process  is  taking  place  in  the  local  order.  In  general  n 
is  expected  to  be  different  in  magnitude  from  p  since  n  is  related  to  the 
local  structural  order,  whereas  p  is  related  to  the  long  range  inter-cluster 
order. 


The  exponential  decay  factor  of  equation  (17)  originates  from  the  cluster- 
dipole  relaxation  process  in  which  all  memory  of  the  initial  prepared  state  is 
lost.  These  are  processes  that  involve  an  external  energy  system  such  as  a 
thermal  bath  and  are  describable  by  a  deterministic  rate  equation  where  &  is 
the  cluster  relaxation  rate.  In  pseudo  one -dimensional  systems  such  as  tfee 
Hollandites  the  cluster  relaxation  rate  is  identical  to  the  recombination 
rate  for  the  effective  cluster  charges,  however  in  two-dimensional  systems, 
such  as  those  of  the  water  adsorbates,  additional  routes  may  be  available  for 
the  cluster  dipole  relaxation  which  do  not  involve  recombination.  The  time 
scale  for  the  relaxation  of  the  cluster  is  then  not  the  same  as  the  time  scale 
for  the  onset  of  the  inter-cluster  exchange  processes. 

3.3  The  composite  process. 

In  the  preceeding  sections  the  presence  of  clustering  in  the  systems 
described  has  forced  the  ion  transport  processes  to  be  divided  into  intra- 
and  inter-cluster  motions  with  distinctly  different  effects.  When  considering 
the  response  of  these  systems  these  two  types  of  motions  have  to  be  combined 
into  a  composite  process. 

Since  we  are  interested  in  the  response  to  a  spatially  uniform  field 
it  is  necessary  to  consider  the  time  development  of  an  initial  (prepared) 
state  with  spatially  uniform  polarisation,  figure  3(c).  The  total 
polarisation  of  this  state  is  identical  to  that  given  by  the  separation  of 
the  effective  cluster  charge,  5,  over  macroscopic  distances  and  these  therefore 
represent  equivalent  descriptions  of  the  total  polarisation  of  the  system  by 
an  applied  field. 


Inter-cluster  exchanges  will  transport  the  effective  charge  through  the 
system  as  described  in  section  3.1  such  that  (nt^)  ^  represents  the  probability 
that  an  effective  charge  is  untransported  at  time  t^  and  { n ( t  -  t.)r  represents 
the  probability  that  an  effective  charge  available  at  time  t^  hasoeen 
transported  at  time  t. 


The  effective  cluster  charge  can  take  part  in  inter-cluster  exchange 
processes  as  a  charged  entity  because  its  separation  from  the  counter  charge, 
£  ,  is  the  cross-over  length  between  motion  as  an  independent  charge  and  as  a 
bSund  cluster  dipole.  The  recombination  of  effective  charges,  cluster 
relaxation,  will  act  as  a  competitive  process  to  the  inter-cluster  exchange 
transport  and  thus  the  probability  of  finding  a  unit  charge  available  for 
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transport  at  time  t^  will  be 

exp{  -  U^)-11  (ntj^)”15 


and  the  contribution  of  an  inter-cluster  exchange  initiated  at  time  t..  to  the 
current  observed  at  time  t  is 


jo  exp{  -  uct1)  Ut^-®  (ntj^)"13  Ct) ( t  -  t^]411 


where  j  is  the  bare  current  element  for  nearest  neighbour  site  ion  hopping 
unperturbed  by  the  internal  structure  of  the  cluster  (Dissado  and  Hill,  1980). 
The  current  observed  at  time  t  is  the  average  of  the  contributions  represented 
by  equation  (20)  ,  and  is  thus 


<J(t)> 


Jo  U/“c) 


r(i+p)  r(i-p) 


C}  f*  e  c  1(w  t1)-nt7p(t-t1)+p  dt.  (21) 
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where  T(l+p)  r(l-p)  t  is  the  normalisation  factor  of  the  distribution, 
figure  6,  and 

i0  (C/»c)'n  -  jc  (22) 

is  the  clothed  (renormalised)  cluster  current.  Carrying  out  the  integral  gives 

<J(t»  -  jc  .  V  «3) 

where  T(  )  is  the  gamma  function  and  -F.(  ;  ;  )  the  confluent  hypergeometric 
function  (Slater,  1960). 

The  complex  dielectric  susceptibility  is  given  by  the  one-sided 
fourier  transform  of  <J(t)>  as 


N  d  r/1  . 

c  c  r(l-p-n) 

c  wckT  (1-in)  r(l-p) 


F (u/u)  ) 
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in  which  d£  («eCc)  is  the  unpertutbed  cluster  dipole  strength,  N  the  number 
density  of  clusters  and  F(u/u  )  the  dielectric  shape  function  asca  function 
of  frequency  scaled  to  the  cluster  relaxation  rate  which  is  given  by 

•  ( rrs] >  (25) 
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in  which  2^(  »  ;  ;  )  is  the  gaussian  hypergeometric  function  (Slater,  1966). 

The  dielectric  function  contained  in  equation  (24)  adequately  describes 
the  response  observed  in  materials  exhibiting  the  anomalous  low  frequency 
dispersion  behaviour  as  can  be  seen  in  figures  7  and  8  where  it  is  applied  to 
results  obtained  for  Hollandite,  humid  sand  and  a  heterogeneous  ionic  glass. 
The  assymptotic  forms  of  F(w/wc)  at  high  and  low  frequencies,  with  respect 
to  u)c  give 


x'(uj)  «  x"(«)  "  X'(»c)  (w/uc) 
X*(w)  «  x"<»>  =  X*  («c)  (w/u.)1 

where  the  factor  x'(“  )  is  given  by 

c 


0)  >  U) 


(26a) 

(26b) 


(27) 


x*(«c) 


j  N 


d 

c 


u  kT 
c 


with  T  the  temperature  and  k  Boltzmann's  constant. 


Kramer s-Kronig  relationships  are  obeyed  throughout  and  in  both  the 
assymptotic  regions  relationship  (2)  is  obeyed,  a  result  which  allows  the 
indexes  p  and  n  to  be  deduced  experimentally  with  some  accuracy  however  closely 
they  approach  either  unity  or  zero.  The  frequency  w  determines  the  cross¬ 
over  from  high  to  low  frequency  behaviour  and  no  sem§lance  of  a  loss  peak  is 
observeable  in  x"(u)  although  a  small  plateau  can  occur  in  x'(w)  for  large 
values  of  p  and  smaller  values  of  n.  This  means  that  a  loss  peak  arising  from 
short  range  hopping  does  not  occur  in  these  materials  although  such  a  loss  peak 
has  been  expected  (Poliak  and  Geballe,  1961),  since  the  site  dipole  hopping 
processes  combine  to  produce  a  whole  cluster  relaxation  which  goes  over 
continuously  to  the  longer  range  inter-cluster  exchange  transport  at  the 
frequency  at  which  the  loss  peak  would  be  expected. 

The  amplitude  of  the  response  {x'(w  ) }  at  the  cross-over  frequency 
o>  simultaneously  represents  a  continuous  cross-over  between  the  dipole 
strength  arising  from  a  cluster  polarisation  of  length  £  and  effective  dipoles 
arising  from  transport  of  the  effective  cluster  charge  over  inter-cluster 
distances . 


For  these  reasons  both  regimes  of  the  anomalous  low  frequency  dispersion 
must  be  regarded  as  a  single  distinct  loss  feature  whose  properties  are  scaled 
to  that  of  a  cluster  such  as  its  length,  £  ,  polarisation  and  recombination/ 
relaxation  rate  o>  .  Only  if  the  cluster  pSssesses  additional  routes  for 
polarisation  relaxation  which  do  not  involve  recombination  will  a  separate 
loss  peak  be  resolvable,  which  nontheless  will  also  form  part  of  the  conq>osite 
system  (Eden  et.al.,  1980) 

3.4  The  relationship  of  anomalous  low  frequency  dispersion  and  the  indices 
p  and  n  to  material  structure. 

The  theoretical  description  in  the  previous  sections  has  shown  that 
anomalous  low  frequency  dispersion  is  a  response  resulting  from  the  relaxation 
of  a  spatially  uniform  displacement  of  charge  in  a  system  with 
aspatially  disordered  partial  occupancy  of  ion  binding  sites.  Initially,  for 
€  <  t  <  u  ,  the  spatially  uniform  displacement  over  macroscopic  distances 

evolves  into^identical  displacements  correlated  over  the  cluster  size  £  .  For 
times  t  >  u  exchanges  of  ions  between  clusters  transports  an  effective  charge 
to  the  boundaries  of  the  system  and  alters  the  internal  structure  of  the 
intervening  clusters  thereby  producing  a  steady  state  distribution  of  their 
structures. 

i)  The  index  n 

This  index  has  previously  been  identified  (Dissado  and  Hill,  1979)  as 
the  dynamic  correlation  index  for  the  intra-cluster  displacement  motions,  it 
therefore  measures  the  efficiency  with  which  displacements  are  transmitted 
within  the  cluster.  An  alternative  definition  of  n  is  obtained  from  the 
average  energy  required  to  establish  the  prepared  state  of  spatially  uniform 
displacement,  as  a  fraction  of  the  maximum  possible  energy  required  to 
establish  this  state,. 

n  ■  E  /c  (28) 

tIC  is  also  equivalent  to  the  maximum  binding  energy  of  the  cluster  when  it  is 
created  from  the  prepared  state  with  E  the  average  binding  energy  (Dissado, 
1982;  Hill  and  Jonscher,  1982).  In  these  terms  small  values  of  n  (  0) 


(12) 


correspond  to  a  weak  distortion  of  the  spatially  uniform  initial  state  and 
hence  a  highly  irregular  cluster  of  distortion  such  as  might  result  from  an 
interstitial  ion.  Large  values  of  n  (-*■  1)  correspond  to  clusters  that  strongly 
distort  the  spatially  uniform  initial  state  and  are  tightly  bound.  Such 
clusters  will  have  a  highly  regular  internal  structure  which  is  nearly 
dissociated  from  the  regularity  of  the  host  lattice  or  matrix  from  which  they 
are  created.  A  value  of  n  identical  to  unity  would  be  impossible  since  it 
would  correspond  to  clusters  which  are  completely  divorced  from  their 
environment. 

The  Hollandites  typically  have  values  of  n  in  the  range  0.6  to  0.8,  which 
shows  that  the  ion-ion  repulsions  which  cause  the  clustering  substantially  over¬ 
ride  the  regularity  of  the  binding  site  array  (Beyeler,  1976).  A  value  of  n 
of  0.8  however  although  high  is  by  no  means  as  high  as  some  that  have  been 
observed  in  loss  peak  processes  (Hill,  1981a;  Shablakh  et.al.,  1982a  &  b)  where 
steric  and  topographical  features  almost  isolate  one  crystalline  domain  from 
another. 

ii)  The  index  p 

This  index  has  been  identified  with  the  dynamic  correlation  index  for 
inter-cluster  exchange  processes  (Dissado  and  Hill,  1979).  Therefore  in  these 
systems  p  indexes  the  efficiency  of  transport  of  the  cluster  effective  charge. 

An  alternative  definition  of  p  also  exists  in  terms  of  the  average  energy 
given  up  by  the  system  when  the  prepared  state  is  converted  to  the  steady 
state  equilibrium  through  inter-cluster  exchange  processes  as  a  fraction  of  the 
maximum  energy  available 

p  ■  E/ri  (29) 

The  probability  density  of  cluster  structures,  figure  6,  shows  that  when  p 
is  small  and  approaches  zero  the  clusters'  structures  remain  almost  identical 
as  in  the  prepared  state  and  the  distribution  is  almost  a  delta  function. 

When  p  is  large  (-►  1)  the  range  of  cluster  structures  is  wide  and  essentially 
unbounded.  Thus  as  p  approaches  unity  the  inter-cluster  exchange  fluctuations 
carry  the  change  of  cluster  structure  far  away  from  their  origin.  In  the 
systems  showing  a  low  frequency  dispersion  the  inter-cluster  exchange  processes 
transport  an  effective  charge,  hence  as  p  approaches  unity  the  effective 
charge  is  transmitted  so  as  to  affect  very  many  clusters.  In  the  limit  of 
p  becoming  unity  the  transported  effective  charge  affects  all  clusters  and 
the  distribution  becomes  unbounded. 

The  same  understanding  of  p  can  be  obtained  if  section  3.1  is  followed 
in  regarding  E  as  the  energy  stored  in  the  structural  fluctuation  defined  as 
the  prepared  state.  This  energy  is  available  through  inter-cluster  exchange 
processes  to  transport  an  effective  charge.  The  maximum  energy  required  by 
inter-cluster  exchanges,  Kn,  is  that  required  to  transport  the  effective 
charge  to  the  sample  surface  and  p  is  the  fraction  of  this  that  is  available 
from  the  structural  fluctuation.  A  value  of  p  identical  to  unity  would 
correspond  to  a  situation  in  which  sufficient  energy  existed  in  structural 
fluctuations  to  completely  remove  the  effective  charge  from  the  sample.  The 
range  of  transport  involved  would  be  infinity  and  result  in  a  normalisation 
factor  of  infinity  in  equation  (16).  Because  an  effective  charge  is  ejected 
this  value  is  impossible  to  realise  unless  the  cluster  reduces  to  a  single 
ion,  in  which  case  a  D.  C.  current  ensues. 

The  values  of  p  typically  observed  in  low  frequency  dispersion  lie 
between  0.7  and  0.99  and  correspond  to  systems  whi^h  approach  perfect  transport. 
The  long  range  structural  description  of  such  systems  is  one  of  a  wide 
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distribution  of  types  of  cluster  structures  with  the  average  clusture 
structure  described  by  the  index  n.  For  example  in  the  Hollandites  we 
may  expect  to  find  a  range  of  differing  clusters  corresponding  to  the  number 
of  ways  of  combining  the  different  supercell  elements  within  the  cluster 
correlation  length.  This  type  of  structure  could  be  expected  to  occur  whenever 
there  exists  a  number  of  local  structural  elements  whose  size  is  incompatible 
with  the  cluster  correlation  length  which  will  be  determined  by  the 
interactions  between  the  elements  arising  from  their  assembly.  This  description 
could  be  expected  to  be  appropriate  to  the  types  of  systems  known  to  exhibit 
the  anomalous  low  frequency  dispersion  response. 

3.5  Hydrogen-bonding  systems 

In  the  preceeding  sections  pseudo  one-dimensional  ionic  systems  such 
as  the  Hollandites  have  been  mainly  used  as  examples  of  low  frequency 
dispersion.  A  two-dimensional  network  of  hydrogen  bonded  molecules  however 
constitutes  a  second  class  of  materials  which  exhibit  this  type  of  response, 
often  with  some  special  features  of  their  own.  In  these  systems  clustering 
will  occur  because  hydrogen-bonding  will  over-ride  the  binding  site  array 
structure  on  a  local  scale.  The  proton  in  the  H-bond  performs  the  function 
of  the  ion  in  the  Hollandites  and  the  cluster  will  be  polarised  by  a 
succession  of  transfers  in  the  H-bonded  system  leaving  a  hydrogen  bonded 
-OH  group  at  the  end  of  one  path  and  a  hydrogen  bonded  -OH_  at  the  other. 

As  with  the  Hollandites  there  will  be  a  number  of  paths  such  that  the  numher 
of  transfers  increases  as  the  individual  charge  separation  energy  decreases. 

In  a  similiar  manner  effective  proton  transfer  over  distances  greater  than 
the  cluster  correlation  length  will  generate  inter-cluster  exchange  processes. 
Because  only  nearest  neighbour  proton  hops  are  allowed,  direct  proton  transfer 
being  inoperative  except  as  a  thermally  activated  mechanism  in  a  completely 
connected  system  (Freund,  1981),  a  two-dimensional  spatiality  is  necessary 
to  give  the  limited  number  of  alternative  routes  required  for  the  power 
law  'evolution'  mechanism. 

The  intra-cluster  regularity  of  these  systems  can  be  expected  to  be 
less  than  that  in  Hollandites  because  of  the  presence  of  weak  as  well  as 
strong  bonds  and  this  is  reflected  in  the  much  lower  value  of  n  that  is  often 
observed,  typically  of  the  order  of  one  half.  Since  the  molecular  structure 
can  be  expected  to  distort  somewhat  to  accomodate  the  incompatibility  of  the 
short  range  H-bonding  and  long  range  site  structure  the  index  p  should  be 
slightly  less  than  those  observed  in  the  Hollandites  and  is  typically  of  the 
order  of  0.8  to  0.9. 

The  main  difference  that  can  be  observed  in  these  systems  is  the 
presence  of  a  subsidiary  loss  peak,  as  shown  in  figure  9,  which  cannot  be 
separated  from  the  onset  of  the  low  frequency  dispersion  even  although  the 
loss  peak  frequency  can  be  shifted  through  several  orders  of  magnitude  by 
variation  of  one  of  the  several  parameters  such  as  temperature  or  humidity. 

This  behaviour  occurs  because  the  two-dimensional  system  has  additional  routes 
for  relaxation  which  do  not  involve  charge  recombination.  In  this  case  the 
characteristic  frequency  for  the  anomalous  dispersion,  w  ,  is  less  than  that 
of  the  loss  peak,  to  ,  and  the  cluster  dipole  has  relaxed  even  although  the 
effective  charge  is^still  available  to  take  part  in  the  inter-cluster  exchanges. 
The  response  can  therefcse  be  regarded  as  the  sum  of  separate  low  frequency 
dispersion  and  loss  peak  processes,  however  in  these  systems  one  further 
feature  leaves  them  coupled,  the  temperature  dependence  of  us  and  u>c  is 
identical.  This  can  be  understood  if  the  recombination  is  regarded  as  a  direct 
intra-cluster  transfer  of  a  proton  by  activation  over  a  H-bond  potential  barrier 


t1 
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with  the  pre-exponential  factor  a  proton  vibration  frequency  multiplied  by 
an  attempt  factor  giving  the  total  number  of  different  routes  available  for 
recombination  (Hill  and  Dissado,  1982).  The  relaxation  of  the  cluster  dipole 
will  also  require  an  identical  activation  energy  but  will  include  routes  by 
which  the  molecules  are  rotated  following  H-bond  breaking  which  is  later 
re-created.  Since  the  rotation  motions  are  coupled  to  the  translation  motions 
in  a  H-bond  (Witkowski  and  Wojcik,  1973)  the  same  activated  state  can  be 
entered  through  either  type  of  motion  and  the  phonon  frequency  will  be 
identical  in  oi  and  to  ,  but  u  will  contain  the  totality  of  routes  whereas  u> 
will  only  contlin  thecfractio£  that  leads  to  charge  recombination,  that  is  c 


Number  of  routes  giving  recombination 
Total  number  of  routes  leading  to  relaxation 


(30) 


Under  these  circumstances  a  cluster  loss  peak  is  observed  at  u  whose 
high  frequency  tail  (u>n  )  contains  the  index  n  appropriate  to  the  regularity 
of  the  site  dipole  structure  which  is  also  the  intra-cluster  regularity  of 
the  H-bond  structure.  On  the  low  frequency  side  of  the  loss  peak  (u  )  the 
inter-cluster  exchange  process  starts  to  produce  a  distribution  of  cluster 
structures  but  does  not  yet  lead  to  an  effective  transport  of  charge  because 
the  charges  have  not  yet  recombined  in  the  polarised  intermediary  clusters 
along  the  transport  path,  figure  3.  In  this  frequency  region  the  inter- 
clusti-  exchange  processes  play  the  role  previously  assigned  to  them  in  loss 
peak  /.ocesses  (Dissado  and  Hill,  1979;  Dissado,  1982).  The  loss  peak  index 
m  will  therefore  have  the  same  value  as  that  of  p  in  the  low  frequency 
dispersion  observed  at  frequencies  less  than  u>  ,  as  seen  in  figure  9.  The 
effect  of  humidity  also  finds  a  simple  explanation  in  these  terms  since 
adsorption  of  extra  molecules  will  occur  in  the  vicinity  of  the  primary 
adsorption  sites  expanding  the  cluster  and  rapidly  increasing  the  number  of 
routes  available  for  relaxation.  The  local  structural  order  is  determined  by 
the  short  range  H-bonding  which  is  unaffected  leaving  n  and  the  ratio  of 
frequencies,  equation  (30),  unchanged. 


The  relationship  of  the  cluster  relaxation  process  to  the  low 
frequency  dispersion  is  illustrated  in  a  sequence  of  results  obtained  from 
the  cyclo-alcohols  (Shablakh  et.al.,  1982b).  In  the  plastic  crystal  phases  of 
cyclopentanol  and  cyclohexanol  hydrogen  bonding  over-rides  the  orientational 
ordering  of  the  molecules  and  the  clusters  are  combined  into  a  three-dimensional 
network.  Because  of  the  three-dimensional  spatiality  a  substantial  number 
of  fully  connected  routes  will  exist  allowing  an  activated  direct  proton 
transport  and  D.C.  characteristics.  The  larger  number  of  cluster  relaxation 
routes  produces  a  separate  relaxation  peak,  figure  10,  whose  temperature 
behaviour  is  the  same  as  that  of  the  D.C.  conductivity.  The  inter-cluster 
exchange  routes  that  are  not  involved  in  the  D.C.  conductivity  act  as  structural 
fluctuations  on  the  loss  peak  giving  an  index  m  »  0.85,  their  low  frequency 
dispersion  contribution  being  short  circuited  by  the  D.C.  current.  Following 
a  time  dependent  transformation  into  an  ordered  crystal  phase  the  three- 
dimensional  network  is  disconnected  and  the  remaining  two-dimensional  boundary 
regions  now  contribute  an  anomalous  dispersion  with  a  value  of  p  identical  to 
that  of  m  for  the  plastic  crystal.  The  values  of  n  in  the  two  states  is  also 
identical  as  can  be  seen  in  figure  11.  A  similiar  behaviour  is  observed  in 
cyclohexanol  where  in  the  ordered  crystal  phase  III  the  cluster  loss  peak  is  not 
quite  resolved  from  the  low  frequency  dispersion  to  which  it  retains  its 
relationship  over  four  decades  of  frequency  change  for  figure  12. 


A.  Limits  of  the  Low  Frequency  Dispersion  Response 


The  mechanism  of  low  frequency  dispersion  described  in  the  previous 
section  relates  the  response  to  the  inter-play  between  short  range  and  long 
range  motions  which  are  scaled  to  the  cluster  size  £  .  This  result  will 
clearly  break  down  at  very  high  frequencies  when  only  the  motions  of  ions 
which  do  not  affect  the  environment  will  respond,  and  also  at  very  low 
frequencies  when  the  ions  may  break  free  from  the  inter-cluster  exchange  routes 
or  the  finite  size  of  the  sample  have  an  effect.  Each  of  these  features 
defines  a  limiting  behaviour  which  may  be  observed  experimentally. 

A.l  High  frequency  limit 

For  frequencies  u  >  u  the  response  of  the  system  is  contributed  by 
the  individual  processes  generated  by  the  perturbations  V  (t)  that  reorganise 
the  structural  ordering  of  the  prepared  states.  The  low  energy  (EE)  motions 
are  those  in  which  individual  site  displacements  occur  in  phase  throughout 
the  cluster.  Their  limiting  intra-cluster  energy  is  equivalent  to  the 
recombination  frequency  u  which  thus  defines  a  cut-off  to  the  infra-red 
divergent  behaviour  as  shown  in  equation  (17).  Higher  frequencies  which 
correspond  to  higher  energy  adjustments  involve  displacements  whose  phase 
can  be  defined  by  a  wavevector  q  that  progressively  increases  (Hohenberg  and 
Halperin,  1977).  At  very  high  frequencies  and  wavevector  only  nearest 
neighbour  motions  are  coupled  together  strongly.  The  effective  charge  will 
therefore  vibrate  as  a  site  bound  ion  in  the  total  cluster  potential.  Its 
vibrations  will  be  equivalent  to  oscillations  of  the  effective  cluster  dipole, 
figure  3(d),  and  their  frequency  A  will  be  a  measure  of  the  binding  energy 
of  the  cluster  dipole  as  a  correlated  entity,  that  is  the  maximum  cluster 
binding  energy.  This  frequency  will  therefore  represent  the  upper  limit 
of  the  infra-red  divergence  and  will  produce  the  oscillations  shown  in  figure 
A  which  become  more  pronounced  as  n  approaches  unity  and  the  average  cluster 
structure  becomes  more  ordered  giving  more  correlated  motion. 


The  high  £r^quency  limit  C  therefore  determines  the  cross-over  from 
the  power  law  oo  response  to  a  different  type  of  behaviour.  If  we  examine 
the  susceptibility  at  this  frequency  we  find  that 

X’(m-C)  »  X'(wc)  (c/«c)n_1 

■  (W  (Nc  dc)  (“c/;) 

■  N  d  /(ckT)  (31) 

o  s  s 

where  N  is  the  number  density  of  site  dipoles  d  . 
s  s 

Expression  (31)  shows  that  c  represents  the  cross-over  from  the  infra-red 
divergent  tail  of  .a  relaxation  process  to  a  resonant  behaviour  with  jo 
representing  the  rate  of  change  of  the  site  dipole  and  therefore  incorporating 
the  damping  factor  (u>  )  of  the  oscillations.  More  precise  details  of  this 
frequency  region  are  Necessary  if  a  complete  description  is  required,  however 
an  echo  of  it  appears  in  the  coupling  strength  (u  )  where  the  coupling 
strength  of  the  high  frequency  process  has  been  cloihed  to  become 

X’(u»c)  -  j0(Nc  dc)  (c/w.)**  («c  kT)_1  (32) 

If  we  can  take 


j 


o 
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as  expected  from  the  perfect  correlation  of  figure  3(d),  then 
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and  the  factor  (u  /£)  defines  the  reduction  of  the  cluster  charge  to  its 
effective  value  by  the  clothing. 

4.2  Low  frequency  limit 

(i)  Sample  size  effect 

For  frequencies  ui  <  u)  the  anomalous  low  frequency  dispersion  (w  p) 
results  from  the  response  ofcthe  I-C  exchange  processes.  As  the  frequency  is 
reduced  lower  energy  inter-cluster  exchange  processes  which  coherently  transport 
a  proportionally  larger  fraction  of  the  effective  cluster  charge,  equation  (8), 
contribute  to  the  response.  A  macroscopic  dipole  moment  is  produced  across 
the  sample  thickness  which  increases  as  the  frequency  is  lowered,  (x'  *  ui  P) . 

A  limit  to  the  polarisation  will  be  reached  for  finite  sized  samples  when  the 
effective  charge  is  completely  transported  to  the  boundaries  thus  defining  the 
lowest  frequency  of  the  I-C  exchange  processes  which  will  correspond_|o  in- 
phase  transport  across  the  sample  thickness  D,  with  wavevector  q  ^  D  . 

The  capacitance  can  be  expected  to  saturate  at  this  frequency,  as  observed  in 
figure  13.  At  lower  frequencies  non-linear  effects  may  occur  (Garton,  1939) 
and  a  D.C.  conductivity  may  arise  if  the  presence  of  an  effective  charge 
allows  the  surface  ions  to  be  discharged. 

(ii)  Onset  of  particle  diffusion 

The  presence  of  an  effective  charge  transport  by  way  of  the  low  frequency 
dispersion  mechanism  relies  on  the  restriction  of  ion  or  proton  motions  to 
their  respective  one  or  two-dimensional  pathways  along  which  structural 
readjustments  can  transport.  When  the  ions  or  water  molecules  can  be  removed 
completely  from  the  pathways  they  will  become  free  to  diffuse  in  the  three- 
dimensional  host  matrix  independently  of  the  clusters  which  are  confined  to 
the  paths.  The  mechanisms  causing  such  scattering  will  need  sufficient  energy 
to  remove  an  ion/molecule  from  its  binding  site  and  the  cluster  in  addition  to 
a  change  in  momentum  vector  and  can  be  composed  of  either  activated  phonon 
processes  or  scattering  off  topographical  imperfections  of  the  host  matrix. 

A  time  t^  can  therefore  be  defined  as  the  lif|time  with  respect  to  these 
processes.  At  frequencies  lower  than  V t,  the  response  assumes  the  familiar 
three-dimensional  diffusion  form  (Carslaw  and  Jaeger,  1959) 

X'(o))  *  x"<u>)  ■  w  ^  ;  u>  <  u>d  (35) 

The  cross-over  frequency  will  define  a  path  length  £,  for  ion/molecule  three- 
dimensional  scattering  and  a  cross-over  of  this  form  nas  been  discussed 
(Shapiro  and  Abrahams,  1981)  for  the  transport  of  electronic  w^gepackets.  For 
the  latter  case  has  been  estimated  to  be  of  the  order  of  10  Hz  or  higher, 
a  frequency  where  it  would  be  difficult  to  observe  the  cross-over.  In  the 
ionic  systems  discussed  here,  however,  it  can  be  expected  that  the  cross-over 
frequency  would  be  much  lower  and  in  many  materials  such  as  the  Hollandites 
with  their  weak  interaction  between  channels  it  may  not  occur  at  all  in  the 
observeable  frequency  range.  It  has  however  been  observed  in  several  systems 
including  ionic  glasses,  pyridine  (Shablakh  et.al.,  1982a)  and  the  organic 

membrane  egg  lecithin  with  cholestrol  (Jonscher,  1978).  The  last  of  these 
examples  is  shown  in  figure  14  where  the  cross-over  is  shown  in  a  striking 
fashion  and  the  diagram  reveals  in  addition  that  the  temperature  behaviour  of 
Uj  is  not  related  to  that  of  u  .  We  may  expect  w.  >  at  high  temperatures 
in  this  system  when  the  clusters  themselves  would-be  dissociated  by  the  free 
particle  diffusion,  perhaps  indicating  the  onset  of  biological  deterioration 
in  this  particular  case. 
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The  onset  of  particle  diffusion  carries  the  system  into  the  classical 
macroscopic  regime  and  although  lower  frequency  processes  may  be  observed 
they  will  not  be  discussed  here. 


5.  The  Relationship  of  Low  Frequency  Dispersion  to  Power  Law  Electrical  Noise 


Expression  (24)  for  x’(u)  can  be  converted  to  an  A.C.  conductivity 
cr( oi)  as,  (Kubo,  1957), 

a(m)  -  <i)XH(u)  (36) 

giving 

a(ai)  s  x'(wc)  «£  -lp  *  x'(“c)  i  “  <  uc  (37a) 

and  o(w)  =  x'^  )  u>n  ;  to  >  <o  (37b) 

c  c  c 


In  materials  exhibiting  low  frequency  dispersion  p  is  often  very  close  to  unity 
and  hence  cr(co)  in  the  region  u  <  u  has  often  been  assumed  to  be  the  D.C. 
conductivity,  from  which  it  is  virtually  indistinguishable  as  shown  in  figure 
15.  However  a  simultaneous  measurement  of  x'(o>),  which  has  the  same  frequency 
dependence  as  x"(m),  will  reveal  that  this  is  not  the  case.  The  relationship 
between  A.C.  conductivity  and  D.C.  conductivity  described  by  Almond  et.al. 
(1982)  for  Na  8-alumina  results  from  equation  (37)  with  p  *  0.99  and  x'(wc) 
essentially  independent  of  <o  .  This  latter  result,  which  is  not  of  general 
applicability  (Jonscher,  197§),  has  implications  for  the  temperature 
dependence  of  c  and  j  in  Na  B-alumina  and  related  materials. 


That  low  frequency  dispersion  does  not  represent  a  D.C.  conductivity 
can  be  seen  from  its  derivation  since  the  effective  cluster  charge  that  is 
transported  decays  via  cluster  relaxation  and  recombination  giving  an  average 
current  <j£t)>  in  equation  (23)  that  clearly  decays  monotonically  with  time 
for  t  >  ui  ,  albeit  extremely  slowly,  figure  16.  The  behaviour  of  <J(t)>  gives 
a  description  of  the  processes  involved  in  low  frequency  dispersion.  Initially 
a  steep  decay  (^t  ;  results  as  the  site  displacement  of  a  unit  charge  is 

clothed  by  the  cluster  readjustments  to  give  the  nett  displacement  of  an 
effective  charge  over  the  cluster  length  £  .  The  effective  charge  will  be 
smaller  as  the  cluster  binding  energy  increases  resulting  in  a  larger  value 
of  n  which  is  necessary  to  achieve  the  larger  reduction  in  current  required 
over  the  cluster  size.  At  times  greater  than  t  =  oo  ,  which  is  an  intra¬ 
cluster  transit  time  for  the  effective  charge,  inter-cluster  exchanges 
transport  the  effective  charge  through  the  disordered  array  of  clusters  with 
an  efficiency  denoted  by  p,  such  that  no  further  decrease  of  charge  occurs 
when  p  is  unity.  The  hypothetical  case  of  n  *  0  and  p  *  1  will  therefore 
correspond  to  a  perfect  D.C.  current. 


When  the  electrodes  to  the  system  have  such  a  nature  that  it  is  possible 
to  exchange  charge  with  the  ionic  system  without  deteriorating  the  electrodes 
or  building  up  a  neutral  barrier  a  D.C.  current  will  flow  under  constant 
operting  conditions  whose  magnitude  will  be  dependent  upon  the  sample  size 
since  the  t^  process  will  be  cut-off  at  a  time  equivalent  to  the  sample 
transit  time  for  the  effective  charge.  The  processes  giving  rise  to  the 
low  frequency  dispersion  will  now  be  observed  as  a  noise  on  the  current.  We 
have  previously  related  expression  (16)  to  power  law  noise  in  electrical 
systems  (Hill  et.al.,  1981)  with  the  result 

Noise  Power  <*  (l/w)a  ;  a  ■  2(l-p)  (38) 

The  present  description  of  low  frequency  dispersion  shows  that  the  noise  arises 
from  structural  order  fluctuations  in  the  disordered  array  of  clusters  by  means 


of  which  an  effective  charge  is  transported  as  a  current.  Energy  is  lost 
irreversibly  when  these  fluctuations  drive  the  thermal  processes  involved 
in  intra-cluster  relaxation  as  can  be  seen  from  the  decay  of  response 
current  <J(t)>  ,  as  shown  in  figure  16.  Thus  intra-cluster  relaxation  also 
contributes  a  thermal  noise  to  the  system. 

The  power  index,  a,  of  the  power  law  noise  is  dependent  upon  p,  such 
that  white  noise  results  if  p  is  unity  and  the  correlation  of  transport 
processes  perfect.  This  corresponds  to  an  unbounded  steady  state  distribution 
of  cluster  structures.  When  p  is  zero  the  cluster  distribution  is  a  delta 
function  corresponding, to  a  superlattice  of  identical  clusters  and  the  noise 
becomes  gaussian  (=  u  ) .  The  distinction  between  these  two  limits  lies  in 
the  number  of  random  inter-cluster  exchange  events  in  a  fixed  time  interval 
(Feynman  and  Hibbs,  1965)  given  by  the  normalisation  factor  of  equation  (16) 
which  is  infinite  for  the  unbounded  distribution  producing  the  white  noise 
and  unity  for  the  superlattice  which  exhibits  gaussian  noise. 

These  limits  are  unachievable  in  practice  for  the  systems  considered 
here  and  values  of  a  in  the  range  from  about  0.4  to  zero  are  to  be  expected. 
When  the  classical  diffusion  process  occurs  we  may  consider  p  to  be  0.5,  in 
which  case  a  is  unity  and  the  noise  becomes  inversely  proportional  to  the 
frequency.  In  electronj.g  system^*11  which  the  diffusion  process  starts  at 
very  short  times  (<  10  to  10  s)  it  will  be  this  process  that  will  be 
commonly  observed,  whereas  in  systems  such  as  polymers  we  may  expect  values 
significantly  different  from  unity  (Pender  ans  Wintle,  1979) . 

It  should  be  noted  that  a  normal  D.C.  response,  equation  (3),  will  be 
obtained  for  systems  possessing  electrodes  that  pass  a  current,  the  anomalous 
low  frequency  response  only  being  observed  when  the  electrodes  block  the 
passage  of  a  current,  or  the  system  is  confined  to  structural  domains  which 
are  not  connected  to  the  sample  electrodes. 

6.  Relationships  to  Other  Models 

The  description  of  the  low  frequency  dispersion  mechanism  presented 
here  relates  the  response  to  the  detail  of  the  material  structural  ordering 
on  a  size  scale  which  links  the  microscopic  motions  of  the  local  structure 
with  the  near  macroscopic  behaviour  of  the  whole  structure.  The  technique 
adopted  to  achieve  this  result  is  in  the  spirit  of  the  path  integral 
formalism  (Feynman  and  Hibbs,  1965).  A  specified  fluctuation  of  the 
structural  order  is  prepared  as  a  wavepacket  of  the  weighted  microscopic 
disordering  process  and  allowed  to  evolve.  In  this  manner  dynamic  indices 
determining  the  evolution  can  be  related  to  the  degree  of  structural  change 
arising  in  re-creating  the  disordered  steady  state  equilibrium  state  from 
the  prepared  state. 

This  method  is  by  its  nature  semi-microscopic  and  semi-stochastic  and 
hence  wide  ranging.  Several  models  however  have  been  proposed  for  this  and 
related  phenomena  which  cover  some  of  the  features  described  here  but  are 

insufficiently  comprehensive  to  describe  them  all.  These  will  be  briefly 
described  so  that  their  relationship  to  the  present  work  can  be  made  explicit 
and  their  differences,  if  any,  outlined. 

6.1  Jonscher's  screened  hopping  model 

In  this  approach  it  was  suggested  (Jonscher,  1975)  that  subsequent  to 
a  hopping  transport  a  charge  carrier  was  screened  by  its  environment  reducing 
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Che  charge  Co  a  facCor  £  of  ics  original  value.  As  a  resulc  elemencary 
cons ideraC ions  o£  scored  and  dissipaced  energy  gives 

cot (mr/2)  -  £r(u)l> '-“eM  "  (1  "  f)/f  (39) 

In  chis  model  che  incra-cluscer  region  of  response  corresponds  Co  Che  motion 
of  weakly  screened  charges  (f  -*■  1) ,  which  score  a  large  fraccion  of  che 
elecCrical  energy  in  a  recoverable  manner.  This  descripcion  is  equivalenc 
Co  chac  of  a  clusCer  whose  binding  energy  and  inCemal  correlacion  of  mocions 
allows  it  Co  score  a  large  proporCion  of  energy  in  Che  relacive  motions  of 
ions.  However  because  Che  ions  are  noC  screened  from  each  ocher  these  motions 
are  strongly  coupled  within  the  cluster,  as  recognised  by  Beyeler  (1976)  and  - 
not  weakly  interacting  as  suggested  by  Jonscher.  The  weakness  of  the 
screening  in  his  model  arises  from  the  poor  binding  to  the  counter-ion 
potential  within  che  cluster. 

If  we  denote  p  as  1  -  n'  Jonscher 's  model  is  again  physically  accurate 
in  describing  the  inter-cluster  exchange  regime  as  arising  from  Che  transport 
of  strongly  screened  charges.  The  weakness  of  the  original  argument  lies  in 
its  divorce  from  the  detail  of  material  structure.  This  was  partly  remedied 
by  Ngai  et.al.  (1979)  who  attributed  the  power  law  response  to  excitations  of 
states  characterising  the  correlations  of  the  particles.  No  attempt  was  made, 
however,  to  relate  this  behaviour  to  a  detailed  structure.  There  is,  in  fact, 
a  considerable  difference  of  interpretation  between  Ngai's  model  and  the 
approach  described  here  and  elsewhere  (Dissado,  1982)  since  in  the  latter  the 
energy  changes  occur  as  alterations  in  the  local  vibrational  frequencies  as 
an  ordered  structure  evolves  to  a  disordered  one  via  correlated  site 
displacements  rather  than  the  excitation  of  a  sub-system  of  states.  An 
additional  feature  is  the  conservation  of  total  energy  between  the  forces 
favouring  ordering  and  those  which  disorder.  Ngai's  suggestion  that  the 
low  frequency  dispersion  involved  two  totally  independent  processes,  the  lower 
frequency  behaviour  dominating  because  of  a  stronger  coupling  constant,  also 
cannot  be  upheld  in  the  light  of  the  experimental  evidence  of  such  features 
as  diffusion.  As  a  consequence  the  relationship  between  the  intra  and  inter¬ 
cluster  regimes  was  not  recognised. 

6.2  Cluster  self -energy 

In  several  problems  concerning  a  region  of  structural  order  the  static 
and  dynamic  response  can  be  related  to  the  self-energy  of  the  structurally 
ordered  cluster,  giving  (Hohenberg  and  Halperin,  1977) 

Cx(<l"0*w)3  1  -  Cxo(q-0j<»»)3  1  +  £(q-0,w)  (40) 

where  x  i-3  the  unclothed  susceptibility  and  £  the  self-energy.  In  problems 
involving  a  non-conserved  property  such  as  is  the  case  following  the 
introduction  of  relaxation  into  the  low  frequency  dispersion  mechanism  the 
self-energy  has  a  logarithmic  divergence 

£(q»0,w)  *  In  id  (41) 

and  self-consistent  solutions  are  sought  by  way  of  a  renormalisation  group 
technique  or  exponentiation  (Hohenberg  et.al.,  1972)  leading  to  a  power  law 
dependence  of  the  clothed  susceptibility 

x'(q*0,w)  «  x"(q-0,u)  «  un  1  (42) 

This  behaviour  is  identical  with  that  of  equation  (26)  which  s  suits  from  a 
method  that  assumes  exponentiation  in  that  the  influence  of  higher  order 
cumulants  in  equation  (S)  are  neglected  and  treats  the  problem  in  the  time 


(20) 


rather  than  in  the  frequency  domain. 

The  self-energy  in  equation  (40)  results  from  a  model  of  the  static 
cluster  response  in  which  the  possible  cluster  displacements  are  averaged 
over  the  cluster  sites  to  give  a  site  probability  density  of  displacement 
(Kogut  and  Wilson,  1972)  which  is  non-gaussian.  Values  of  n  approaching 
unity  as  appropriate  to  the  intra-cluster  structure  are  equivalent  to  a  well 
ordered,  highly  correlated,  structure.  When  p  approaches  unity  (n'  ■+  0)  the 
structural  order  is  very  weak  corresponding  to  the  inter-cluster  ordering. 

It  is  interesting  that  the  often  used  modulus  representation  of  dielectric 
permittivity  (cf.  Wong  and  Angell,  1976)  presents  the  experimental  results  in 
the  form  of  equation  (40).  In  the  low  frequency  region,  it  is  shown  in  the 
modulus  plot  of  figure  17,  that 

Real  Cx(q*0,d))]  *  -*■  0  (.43) 

as  u)  -*•  0  and  thus  the  effective  charges  in  the  intra-cluster  exchange  regime 
are  highly  screened  from  each  other  as  we  have  determined.  The  modulus,  or 
inverse  complex  susceptibility,  has  the  interesting  feature  that  the  plot 
of  low  frequency  dispersion  data  yields  a  peak  in  the  imaginary  component 
making  determination  of  relatively  easy. 

6.3  Scaling  approach 

The  scaling  technique  has  been  used  extensively  to  describe  systems  in 
which  order  exists  only  over  a  finite  sub-macroscopic  size.  Basically  the 
aim  is  to  determine  the  extensive  property  (P(L})  of  a  system  as  a  function  of 
its  size.  Starting  from  a  small  sized  element  for  which  a  nearly  exact 
calculation  can  be  made,  the  size  is  increased  by  a  constant  factor  b  and  the 
property  P{bL}  recalculated.  When  the  smaller  system  has  the  same  functional 
dependence  on  size  as  the  larger  the  original  element  can  be  seen  to  be  a 
smaller  sample  of  the  larger  system  with  the  same  intensive  properties.  The 
fundamental  hypothesis  of  the  technique,  which  was  suggested  by  Kadanoff 
(1966),  is  that  a  minimum  element  of  structure  exists  whose  form  is  retained 
as  the  size  scale  increases.  This  assumption  leads  to  the  relationship 

P{bL)  -  f(b,P{L})  (44) 

which  has  been  applied  by  Abrahams  et.al.  (1980)  to  the  analysis  of  the  static 
electronic  conductivity  of  disordered  systems.  The  size  dependence  of  the 
required  property  can  be  determined  from  equation  (44)  by  dimensional 
arguments  only  if  the  system  is  perfectly  ordered,  any  disorder  contributes 
an  infra-red  divergent  self-energy  which  leads  to  an  anomalous  dimensionality 
(Wilson  and  Kogut,  1974;  Abrahams  et.al.,  1980).  This  index  determines  the 
structural  relationship  of  the  scaled-up  element  to  its  precursor  and  is 
directly  related  to  the  dynamic  index  n  (Hohenberg  and  Halperin,  1977)  which 
therefore  describes  the  structural  regularity  in  this  sense.  It  should  be 
regarded  as  a  necessary  factor  in  the  description  of  structural  ordering  in 
addition  to  a  correlation  length. 

When  the  method  is  applied  to  the  present  problem  the  static  susceptibility 
is  first  determined  as  a  function  of  the  length  L.  Determination  of  the 
eigenfrequency  w(L)  for  a  fluctuation  of  the  same  length  then  enables  the 
dynamic  susceptibility  to  be  obtained  (Shapiro  and  Abrahams,  1981)  as  a 
function  of  frequency.  The  field  is  considered  to  couple  individually  to  the 
fluctuation  with  the  static  susceptibility  for  a  given  length  over  a  continuum 
of  lengths  up  to  the  sample  size  D.  If,  for  example. 
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x(L)  ^  LX  and  (u(L)  ^  (45) 

then  x(“)  ^  <un  ^  with  n  ■  1  -  x/y  (46) 

Previous  work  has  been  concerned  with  the  transport  of  an  electron 
wave  packet  (Thouless,  1973;  Shapiro  and  Abrahams,  1981)  and  not  with  wave 
packet  formation  which  would  be  equivalent  to  the  intra-cluster  processes 
considered  here.  The  cluster,  however,  can  be  scaled  in  the  same  way  with 
the  minimum  structural  element  being  the  site  dipole.  The  observation  of  a 
value  of  n  approaching  unity  for  a  frequency  range  of  several  decades  shows 
that  the  susceptibility,  and  thus  the  polarisation  of  the  fluctuation,  is 
only  weakly  increasing  with  fluctuation  size  within  the  cluster,  while  the 
frequency  is  strongly  size  dependent.  The  polarisation  coupling  to  a 
uniform  field  is  that  of  an  in-phase  displacement  and  will  be  independent  of 
the  cluster  size  if  the  site  displacements  are  fully  correlated  as 
independent  fourier  components.  When  the  site  displacements  have  zero 
correlation  the  polarisation  will  be  proportional  to  the  cluster  size.  All 
real  systems  will  lie  between  these  two  extremes  with  n  approaching  unity 
corresponding  to  strong  correlation  as  concluded  in  section  3. 

The  wave-packet  transport  behaviour  is  equivalent  to  the  inter-cluster 
exchange  regime,  and  the  value  of  p  -*■  1  observed  over  a  large  frequency  range 
indicates  a  strong  length  dependent  increase  of  the  polarisation  and  hence  a 
nearly  uncorrelated  inter-cluster  structure.  It  follows  that  the  current 
appropriate  to  this  regime  will  consequently  show  a  weak  decrease  with  length 
giving  a  conductivity  dependent  upon  the  sample  size  as  long  as  the  low 
spatial  dimensionality  (d  $  2)  is  retained  (Abrahams  et.al.,  1980;  Shapiro 
and  Abrahams,  1981).  The  description  presented  here  is  in  agreement  with 
these  conclusions,  cf  figure  16,  however  the  logarithmic  divergences  in  the 
free  energy  are  exponentiated  to  give  power  law  dependencies  in  time  equivalent 
to  an  anomalous  dimensionality  in  the  spatial  property  (Hohenberg  and  Halperin, 
1977) .  In  this  way  the  logarithmic  frequency  and  length  dependence  of 
Abrahams  et.al.  (1980)  will  be  modified  into  a  power  law  behaviour  (Mott  and 
Kaveh,  1981) . 

6.4  Distribution  of  hopping  times 

Bemasconyi  and  his  co-workers  (see  Beraasconyi  et.al.,  1979)  have 
developed  this  approach  in  a  series  of  papers  in  which  the  Mori  formalism 
has  been  utilised  together  with  a  static  distribution  of  inter-site  hopping 
rates.  It  was  shown  that  when  this  distribution  diverged  for  small  rates 
anomalous  conductivity  ensued  of  the  form  of  equation  (26) ,  with  a  non¬ 
uni  versal  index  p  whose  value  depended  on  the  form  of  the  divergence  in  the 
hopping  rate  distribution.  Alternatively  a  D.C.  conductivity  was  found  when  the 
distribution  was  bounded. 

This  approach  bears  some  superficial  resemblance  to  that  presented  here 
in  that  a  self-consistent  solution  to  equation  (4)  results  in  the  Mori 
formalism  and  the  distribution  used  were  of  the  type  to  produce  an  infra-red 
divergence.  However  their  concept  of  a  static  spatial  distribution  of 
hopping  rates  is  quite  distinct  from  that  of  the  present  work.  The  difficulty 
with  such  a  concept  is  that  it  has  not  been  related  to  the  material  structure 
except  insofar  as  a  Poisson  distribution  of  site  separations  leads  to  a 
similiar  result  in  spin  systems  (Theodoru  and  Cohen,  1979).  Because  of  the 
lack  of  a  relationship  with  structure  the  two  regimes  given  by  equations 
26  (a)  and  (b)  were  not  identified  theoretically.  As  a  consequence  the  low 
frequency  dispersion  regime  in  Hollandites  was  taken  to  be  a  D.C.  conductivity 
(Bernasconyi  et.al.,  1979)  and  only  the  high  frequency  regime  was  analysed. 


Because  the  hopping  rate  distribution  was  assumed  to  be  due  to  a  distribution 
of  activation  energies  a  temperature  dependent  value  of  n  was  obtained  which 
was  predicted  to  become  zero  at  T  *  410  K,  at  which  temperature  the  system 
became  D.C.  conducting  through  a  mobility  transition.  The  experimental  data 
for  the  Hollandites,  figure  7,  show  no  such  temperature  dependence  with  the 
same  low  frequency  dispersion  characteristics  observeable  up  to  863  K. 

The  experimental  data  given  by  these  ^uthorg  in  support  of  their  model 
was  obtained  in  a  small  frequency  window,  10  -  10  Hz,  and  over  a  narrow 

range  of  temperatures,  150  K  to  213  K.  Since  u>  moves  to  higher  frequencies 
with  increasing  temperature  the  low  frequency  dispersion  regime  of  equation 
(26a)  will  move  towards  the  lower  part  of  the  frequency  window  effectively 
reducing  the  observed  value  of  n  and  giving  the  form  of  result  observed. 

If  this  regime  were  taken  to  be  a  D.C.  conductivity  then  a  continuous  cross¬ 
over  would  occur  as  it  moved  into  the  frequency  window.  The  conclusions 
drawn  by  these  authors  with  respect  to  the  temperature  dependence  of  the 
conductivity  are  therefore  erroneous  and  caused  by  the  limited  range  of 
measurements  made. 

7.  Conclusions 

A  new  model  of  anomalous  low  frequency  dispersion  has  been  presented 
which  relates  the  mechanism  directly  to  the  disordered  material  structure 
which  is  described  in  terms  which  link  the  local  order  and  the  long  range 
order.  The  two  power  law  regimes  in  frequency  have  been  identified  with 
fluctuations  in  the  structural  order  on  a  cluster  scale  and  on  an  inter¬ 
cluster  macroscopic  scale. 

The  behaviour  can  also  be  expected  in  low  dimensional  disordered 
materials  which  possess  a  potential  electronic  conductivity.  In  this  case 
the  cluster  can  be  identified  with  an  electronic  wave-packet  spread  over  a 
number  of  donor  sites  and  possessing  an  internal  structure  due  to  the  local 
order.  The  polarised  cluster  can  therefore  be  thought  of  as  an  intermediate, 
rather  than  a  Wannier,  exciton.  Since  the  anomajLgus  transport  will  be  cut¬ 
off  at  long  times  by  ele^ron  diffusion  (t  >  10  s)  and  at  short  times  by 
quantum  effects  (t  <  10  s)  the  restricted  time/ frequency  range  available 
will  render  it  likely  to  be  observed  in  electronic  systems  of  low  dimensionality. 
The  long  time  limit  in  such  systems  has  been  shown  to  be  of  the  conventional 
diffusion  limited  form.  The  relation  between  electrical  noise  and  the 
anomalous  low  frequency  dispersion  has  been  examined  and  in  the  diffusion 
limited  case  gives  a  power  noise  rigorously  proportional  to  the  inverse 
frequency.  Variations  from  this  characteristic  can  be  taken  as  evidence  for 
a  non-perfect  diffusion  in  the  system  un^er  examination. 

It  has  been  clearly  demonstrated  for  non-electronic  charge  carriers, 
such  as  protons,  that  the  time  scale  over  which  the  anomalous  properties 
will  be  observed  occurs  over  much  longer  times  than  for  the  electronic  case 
and  can  conveniently  be  determined  through  the  noise  power/frequency 
relationship. 

It  has  been  shown  that  even  with  transmitting  electrodes  the  systems 
investigated  here  exhibit  a  sample  size  dependent  conductivity,  decreasing 
with  increasing  size.  The  quantitative  variation  of  the  conductivity  with 
size  has  not  been  investigated  here  as  it  will  depend  on  the  properties  of 
the  specific  system  examined.  It  has  also  been  shown  that  measurement  of 
the  conductivity  alone,  or  the  loss  alone,  may  not  be  sufficient  to 
discriminate  between  low  frequency  dispersion  and  D.C.  conductivity.  It  is 
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considered  that  this  observation  may  have  serious  consequences  in  the 
evaluation  of  solid  fast  ion  conductors  as  power  cells.  The  relation 
between  the  anomalous  dispersion  and  electrical  noise  has  been  examined  and 
found  to  have  a  common  origin  and  related  power  indices. 

Low  frequency  dispersion  is  known  to  occur  in  bio-systems  such  as 
hydrated  proteins  (Bone  et.al.,  1981)  and  the  in-vivo  response  of  plants  to 
electrical  stimulus  (Hart,  1981).  It  can  be  expected  that  the  interpretation 
of  this  phenomenon  which  has  been  presented  here  will  assist  in  a  more 
comprehensive  understanding  of  some  of  the  features  of  these  fields. 

It  requires  to  be  emphasised  that  the  model  developed  here  need  not 
be  restricted  to  electrical  transport.  When  the  inter-cluster  exchange 
processes  carry  no  charge  a  relaxation  loss  peak  will  be  observed 
dielectrically,  but  if  their  motions  are  conjugate  to  an  applied  field  such 
as  pressure  then  a  flow  will  result.  The  separation  of  the  site  motions  into 
intra-cluster  and  inter-cluster  has  been  recognised  by  Cotterill  and  Tallon 
(1980)  who  have  pointed  out  that  the  motion  of  an  atom  along  a  single 
coordinate  in  a  liquid  can  be  simultaneously  a  vibration  and  a  transportation 
We  expect  the  basic  concepts  of  the  model  to  be  of  general  applicability  and 
of  use  in  a  number  of  branches  of  science  dealing  with  imperfect  transport  in 
disordered  systems.  Because  the  model  connects  the  microscopic  structure 
to  the  macroscopic  stochastic  system  with  complete  generality  of  form  it  is 
equally  applicable  to  solids  and  liquids  alike. 
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Figure  1.  Representation  of  anomalous  low  frequency  dispersion  with 
p  *  0.95  and  n  »  0.7.  The  frequency  u  indicates  the  relaxation  rate  and 
the  parallelism  of  the  real  and  imaginary  components  in  this  log/log 
presentation  can  be  clearly  seen  for  both  u  >  u  and  w  <  u  .  The  dotted 
curves  show  the  D.C.  case  and  the  chain  dotted  ?ine  define!  a  frequency 
above  which  the  response  may  be  confused  with  D.C.  behaviour. 

Figure  2.  (a)  Projection  of  the  Hollandite  structure  onto  the  (001)  plane. 

The  barred  ions  are  situated  at  z  *  0.5  and  the  0  and  K  symbols  represent 
the  oxygen  and  potassium  sites,  respectively,  with  the  M  sites  containing 
either  Ti  or  Mg.  (After  Bystrom  and  Bystrom,  1950) 

(b)  The  developed  Hollandite  structure  showing  two  unit  cells  in 
the  z  direction.  One  of  the  channels  within  which  the  potassium  ions  can 
move  in  the  presence  of  ion  vacancies  is  indicated. 

(c)  Representation  of  the  channel  sequences  combined  in  groups 
to  give  the  channel  correlation  length  of  35  A  (Beyeler,  1976).  The  empty 
cells  indicate  vacancies. 

Figure  3.  (a)  Schematic  representation  of  three  inter-cluster  exchange 

mechanisms  which  transport  a  charge  over  the  same  range.  The  arrows  show  the 
ion  transport  and  the  signs  give  the  charge  densities  produced. 

(b)  A  schematic  representation  of  a  spatially  uniform  polarisation 
(ii)  in  a  typical  cluster  (i)  which  indicates  the  equivalence  to  a  whole 
cluster  dipole  with  an  effective  charge  5.  The  ions  and  lattice  counter-charges 
are  shown  for  a  system  with  75%  occupancy.  The  dashed  lines  indicate  the  site 
centres . 

(c)  Schematic  representation  of  all  the  system  processes. 

(i)  Spatially  uniform  cluster  polarisation  ;  (ii)  Some  inter-cluster  exchange 
processes  and  (iii)  Cluster  relaxation  processes. 

(d)  Schematic  representation  of  fully  coupled,  out  of  phase,  ion 
displacements  showing  the  equivalence  to  the  vibration  of  the  effective  cluster 
dipole. 

Figure  4.  The  function  <Kt)  describing  the  time  development  of  the  prepared 
state  for  a  range  of  values  of  p.  As  p  approaches  unity  the  oscillations 
increase  in  amplitude  and  for  p  i  0.9  the  first  oscillation,  at  least,  goes 
negative. 

Figure  5.  A  schematic  diagram  showing  the  decay  and  growth  of  the 
fluctuations.  Initially  it  has  been  assumed  that  the  individual  elements 
are  all  at  their  average  value  in  the  prepared  state  and  the  fluctuations 
give  rise  to  the  steady-state  distribution  shown  on  the  right-hand  side.  The 
diagram  has  been  constructed  for  p  =  0.3  and  indicates  both  the  gaussian 
approach  to  zero  time  and  the  resonance-like  behaviour  for  £t  being  of  the 
order  of  unity  as  shown  in  figure  4. 

Figure  6.  The  probability  density  functions  for  a  range  of  values  of  the 
index  p.  The  fluctuation  scale  has  been  normalised  to  the  average  value  for 
each  distribution. 


Figure  7.  Master  curves  of  susceptibility  as  a  function  of  frequency  for 
(a)  Hollandite  of  composition  I^Al-TigO-j,  (b)  Hollandite  of  composition 
K.  gMg8  gTiy  and  (c)  humid^sand.  The  continuous  functions  have  been 

obtained  from  equation  (25)  with  the  values  of  p  and  n  determined  as 
(a)  p  *  0.92  ,  n  «  0.7  (b)  p  *  0.92  ,  n  =  0.6  and  (c)  p  *  0.8  ,  n  *  0.94. 

The  curves  have  been  displaced  for  clarity.  The  data  for  (a)  and  (b)  was 
taken  from  Jonscher  et.al.,  (1981)  and  for  (c)  from  Shahidi  et.al.  (1975) 

Figure  8.  Master  curve  of  the  susceptibility  of  a  glass  (Allied  base  glass) 
from  the  data  of  Doyle  (1981).  The  continuous  curves  have  been  determined  from 
equation  (25)  with  p  *  0.995  and  n  ■  0.5.  The  diagram  is  scaled  at  800  K. 

Figure  9.  Master  curve  for  the  loss  component  of  susceptibility  for  hydrated 
bovine  serum  albumin  from  the  data  of  Eden  et.al.  (1980).  The  normalisation 
has  been  carried  out  in  terms  of  the  percentage  water  adsorbed  and  the 
diagram  is  scaled  at  8.7  wtZ.  The  normalisation  of  both  the  anomalous 
dispersion  and  the  loss  peak  clearly  indicates  a  common  basis  for  both  processes. 

Figure  10.  The  real  and  imaginary  components  of  the  permittivity  of  the 
plastic  crystal  state  of  cyclopentanol  showing  a  loss  peak  at  high,  relative, 
frequencies  and  pure  D.C.  conductivity  with  no  dispersion  in  the  capacitance 
at  low  frequencies.  The  diagram  is  scaled  at  180  K  and  the  data  is  taken  from 
Shablakh  et.al.  (1982b) 

Figure  11.  The  loss  component  of  the  susceptibility  for  the  crystal  state 
II  of  cyclopentanol  showing  a  weak  loss  peak  at  high,  relative,  frequencies 
which  does  not  normalise  with  the  anomalous  dispersion.  The  plot  is  scaled 
at  110  K  and  the  data  is  taken  from  Shablakh  et.al.  (1982b). 

Figure  12.  The  loss  component  of  the  susceptibility  in  the  crystal  III 

structure  of  cyclohexanol.  Over  the  normalisation  temperature  range  of 

190  K  to  218  K  the  weak  loss  process  embedded  in  the  low  frequency  dispersion 

remains  in  its  relative  position  showing  that  the  two  processes  have  a 

common  basis.  The  data  is  taken  from  Shablakh  et.al.  (1982b).  Scaled  at  218  K. 

Figure  13.  The  permittivity  of  a  doped  glass  (J  22)  in  the  temperature  range 
between  490  K  and  1066  K.  The  normalised  plot  shows  the  lower  frequency 
range  of  an  anomalous  low  frequency  dispersion  with  evidence  of  saturation  in 
both  the  real  and  imaginary  components  at  the  lowest  reduced  frequencies.  The 
continuous  curves  were  obtained  from  equation  (25)  with  p  *  0.92  and  n  *  0.4. 

The  plot  is  scaled  at  1066  K.  (Doyle,  1981) 

Figure  14.  The  permittivity  of  egg  lecithin  measured  at  two  temperatures. 

At  the  higher  temperature  the  onset  of  the  diffusion  controlled  process, 
with  p  ■  0.5,  can  be  seen  together  with  the  distortion  this  introduces  into 
the  real  component  in  the  anomalous  low  frequency  dispersion  region.  The 
continuous  curves  have  been  determined  from  equation  (25)  with  p  *  0.86. 

The  two  sets  of  data  have  been  displaced  for  clarity  (after  Jonscher,  1978). 

Figure  15.  The  anomalous  low  frequency  dispersion  relationship  of  equation 
(25)  plotted  in  the  form  of  A.C.  conductivity  from  equation  (36),  for  a  range 
of  values  of  p.  As  p  approaches  unity  the  conductivity  becomes  indistinguishable 
from  the  frequency  independent  D.C.  form. 


Figure  16.  Schematic  representation  of^the  average  response  current  as  a 
function  of  time  in  the  region  of  t  *  u  .  An  exponential  relaxation  is 
shown  for  comparison  as  is  the  region  of  response  that  leads  to  the  observation 
of  loss  peaks.  The  perfect  inter-cluster  transport  case,  with  p  -  1,  is  the 
0.  C.  like  behaviour  shorn  in  figure  15. 


Figure  17.  Typical  anomalous  low  frequency  dispersion  presented  in  the 
form  of  the  dielectric  modulus  where  M(iu)  ■  {e(to)}  .  The  diagram  illustrates 

the  approach  of  both  the  real  and  imaginary  components  to  zero  at  zero 
frequency  and  the  presence  of  a  peak  in  the  imaginary  component  of  the  modulus. 
This  diagram  should  be  compared  with  figure  1,  which  presents  the  same  information 
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Abstract 

The  experimental  data  from  a  dielectric  investigation  of  a  series 
of  doped  ceramic  ferroelectrics  has  been  examined  in  order  to  determine 
whether  the  data  showed  effects  due  to  the  finite  size  of  the  ordered 
grains.  It  has  been  found  that  the  order  of  the  ferroelectric-paraelectri 
transition  is  affected  by  the  grain  size  and  it  is  pointed  out  that  this 
has  consequences  on  the  mode  structure  of  the  materials.  No  direct 
evidence  was  found  for  correlation  length  limitations  but  the  correlation 
lengths  in  the  samples  examined  was  determined  as  of  the  order  of  thirty 
lattice  spacings. 


Introduction 


The  aim  of  this  study  was  to  examine  experimental  data  in  order  to 
determine  whether  the  effect  of  finite  sample  size  could  be  discerned. 

To  this  end  we  have  chosen  to  look  at  the  available  data^on  a  doped  ceramic 


ano 


ferroelectric  which  had  been  examined  in  this  laboratory  through  i^s 
dielectric  susceptibility,  x(w)»  in  the  frequency  range  between  10  Hz 
10  Hz  over  a  wide  range  of  temperatures  about  the  critical  Curie 
temperature,  T  .  The  reason  for  this  choice  was  two-fold.  Firstly  a  body 
of  theory  exists  which  predicts  that  those  properties  that  formally  diverge 
at  the  critical  temperature,  such  as  the  dielectric  susceptibility  and  the 
relaxation  time,  (w  )  ,  saturate  when  the  correlation  length,  Z,  becomes 

equal  to  the  finite^size,  L,  of  the  sample.  Secondly  the  ceramic  was 
composed  of  finite  sized  grains  with  an  average  diameter  that  was  dependent 
on  the  method  of  preparation  and  of  the  concentration  of  the  uranium  dopant 
that  had  been  introduced.  The  range  of  grain  diameters  available  was  from 
2.5um  to  24pm.  Thus  the  system  afforded  the  possibility  of  examining  the 
predictions  of  the  correlation  length  theory.  The  characteristics  of  the 
samples  used  are  listed  in  Table  I. 


The  theoretical  predictions  can  be  summarised  by  stating  that  when 

Z  >  L  the  correlation  length  in  the  diverging  expressions  can  be  replaced  b 
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where  F(  )  is  a  complex  function  of  the  scaled  frequency  oi/w  .  Expressions 
(1)  to  (5)  give  a  series  of  relationships  that  can  be  examinid  rxperimentally 


Analysis  of  Experimental  Data 

The  experimental  data  showed  tljie  presence  of  two  loss  peaks  in  the 
frequency  window  used  for  observation  .  Superimposed  on  the  loss  peaks 
was  a  substantial  D.C.  conductivity  which,  in  the  frequency  domain,  gives  a 
loss  decreasing  with  increasing  frequency.  The  small  value  of  the 
susceptibility  increment  associated  with  the  loss  process,  Ax'Oo),  and  the 
sensitivity  of  the  magnitude  of  the  increment  to  dopant  concentration 
identifies  the  observed  response  as  arising  from  regions  of  the  samples 
which  are  affected  by  the  impurities.  The  true  ferroelectric  response 
occured  at  frequencies  higher  than  those  measured  but  a  good  estimate  of  the 
ferroelectric  susceptibility  increment  was  obtained  from  the  magnitude  of 
the  dielectric  permittivity  at  10  Hz.  These  values  are  plotted  in  figure  1, 


as  a  function  of  temperature,  for  all  five  samples  that  had  been  investigated 
In  all  but  a  single  case  the  transition  was  found  to  be  of  first  order,  a; 
can  be  seen  by  the  discontinuities  in  the  curves.  The  exception,  sample  J?M35 
was  the  sample  which  had  been  prepared  with  only  a  single  calcining  to  give 
particularly  small  particle  sizes.  The  discontinuity  in  the  traces  for  the 
other  samples  occured  at  th^  critical  temperature.  Making  recourse  to  the 
phenomenological  expression  for  the  free  energy 


G  = 


°io  * ! p2 


b  4  c  6 
+  4P  +  6P 


(6) 


in  which  b  and  c  are  essentially  constant  and  a  is  given  by 


a  '=  4ir(T  -  T  )/C 
o 


(7) 


When  b  is  negative  the  transition  is  first  order  and  the  critical  temperature 
is  given  by 


T 

c 


T  + 
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3Cb 

64rrc 


(8) 


and  there  is  the  discontinuity  in  the  susceptibility  such  that 


X(Tc+)  =  16c/ (3b2) 

and  X(TC-)  =  4c/ (3b2) 

hence  {  x(Tc+) )/ {x(Tc”) }  =  4 


(9) 

(10) 
(ID 


Approximately  this  value  has  been  observed  in  all  the  samples  except 
those  containing  0.4%  uranium,  RM35A  for  which  there  was  no  discontinuity 
and  RM35B  for  which  the  discontinuity  was  smaller,  approximately  two. 


The  case  in  which  no  discontinuous  transition  can  be  seen  can  be 
characterised  by  a  positive  value  of  b  in  equation  (6) ,  and  hence  is  of 
second  order.  In  this  case  there  is  an  infinity  in  x(^®0)  at  T  ,  which 
equal  to  T  ,  of  the  form 


(xO^o)}'1  = 
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H 
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T  >  T 
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(12a) 
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is 


and  the  ratio  of  the  gradients,  d(x)  VdT,  above  and  below  T  has  a  value  of 
two  as  observed  for  sample  RM35A. 


Perovskites,  such  as  Jjhose  investigated  here,  are  expected  to  show  a 
displacive  phase  transition  with  a  soft  mode.  Tne  modification  to  a  first 
order  transition  may  be  attributed  to  electro-stricticn  *  which  arises  from 
the  grain  boundaries  of  the  undamped  crystals  and  leads  to  the  negative 
value  for  b.  It  therefore  appears  that  the  snail  grain  size  in  sample  SM33A 


t 
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alleviates  this  effect  and  allows  the  second  order  transition  to  take  place. 

It  is  generally  considered  that  increases  of  dopant  concentration  decreases 
T  through  a  reduction  of  T  caused  by  lattice  strain.  However  the  dip  in 
T  for  the  larger  partice  sized,  double  calcined,  O.UZ  uranium  sample  (Id-1353) 
might  be  due  to  a  reduction  in  the  magnitude  of  b,  which  then  becomes  positive 
for  the  small  grain  sized,  single  calcined,  sample  of  the  same  concentration 
of  dopant. 


As  stated  earlier  at  lower  frequencies  than  the  ferroelectric  process 
two  loss  peaks  could  be  delineated  wighin  a  background  of  some  D.C.  conductivity. 
An  anomalous  low  frequency  dispersion  which  is  due  to  restricted  mobility  of 
charge  centres  in  the  samples  was  also  observed.  The  two  loss  peaks  were  even 
observed  in  the  unhoped  sample.  It  is  known,  however,  from  the  method  of 
sample  fabrication-  that  the  undoped  material  contains  vacancy  centres  in  the 
perovskite  structure.  Of  the  two  loss  processes  it  was  only  possible  co 
study  the  lower  frequency  one  in  detail.  It  was  found  that  this  process  was 
very  sensitive  to  the  dopant  concentration,  as  indicated  in  figure  2.  It 
seems  likely,  therefore,  that  this  peak  can  be  attributed  to  uranium  atoms 
in  the  doped  materials  and  that  these  lie  at  the  same  type  of  lattice  sites 
as  the  vacancies  in  the  un-aoped  sample.  As  the  low  frequency  loss  peak 
exhibited  critical  behaviour  it  would  appear  that  these  sites  must  lie  withi^ 
the  order.ed  grains  and  are  coupled  to  the  ordering  field  of  the  host  lattice  . 
Insufficient  data  exists  to  allow  any  assignment  to  be  made  to  the  slightly 
higher  frequency  loss  process. 


Scaling  Effects 

’•  i  t 

The  mean  field  theory  outlined  in  equations  (6)  to  (11)  often  gives 
an  approach  to  a  critical  point  such  as  T  with  incorrect  values  for  the 
power  exponents  v,  n  and  y  (*  v{2-q}^.  Different  values  for  these  exponents 
are  predicted  by  more  exact  theories  and  are  generally  called  scaling  effects. 


i  Static  Scaling 

These  are  effects  observed  j.n  the  static  susceptibility,  which  we  have 
taken  to  be  the  permittivity  at  10  Hz  as  regards  the  very  high  frequency 
ferroelectric  process.  The  susceptibility  exponents  for  this  process  are  given 
in  Table  II.  It  can  be  seen  that,  apart  from  the  second  order  transition,  the 
exponents  do  not  have  the  mean  field  values.  Above  T  the  exponent  y  has  a 
value  typical  of  a  three  dimensional  Ising  model,  rather  than  that  of  unity, 
c.f.,  equation  (12).  Below  T  the  value  of  y  in  the  first  order  systems  is 
0.73  independent  of  grain  size  or  of  uranium  concentration.  This  confirms 
our  proposal  that  the  transition  is  first  order  and  of  the  same  dimensionality 
in  all  the  samples  except  those  with  small  grains. 


ii  Dynamic  Scaling  in  the  low  frequency  response 

This  form  of  scaling  establishes  a  relationship  between  the  relaxation 
rate  and  the  susceptibility  increment3  in  the  critical  region.  It  can  be 
succinctly  defined  in  the  following  equations  which  apply  in  the  temperature 
region  in  the  neighbourhood  of  Tc» 


X(w)  ” 


X'('u=0)  F(u)/w  ) 
P 


X'(w)  -  ix"(u) 


(13a) 


x'(u>)  «  x’C^*0)  (w/u  )  (13b « 

X"(“)  *  x'(u=0)  (a)/Wp)n  1  (13c) 

for  u/(i)  >1  and  0  <  n  <1,  with 

P 

X 1  (d>*0)  *  (“p)n  *  and  x'(w=0)  .  *  jT  -.T^  Y  (14) 
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One  further  relationship  has  been  predicted  *  and  should  be  obeyed 
outside  of  the  critical  region,  that  is  at  temperatures  not  close  to  T  . 
This  is  that 

x'Coi-O)  *  (*  )n  ;  with  o>  «  e  A^kT  (15) 

P  P 

The  existance  of  both  relationships  (14)  and  (15)  has  been  revealed^- 
in  the  low  frequency  r^ponse  of  the  samples  investigated  by  the  use  of  a 
diagrammatic  technique  which  constructs  a  master  curve  of  the  response 
with  the  functional  form  of  equation  (13a) .  These  relationships  were 
observed  in  ajl  the  samples,  including  sample  RM35A.  A  typical  locus  of 
1°S t| X C*o=0>  |  }  as  a  function  of  log{(i>  )  J‘}  is  shown  in  figure  3.  The 

inverse  values  in  this  plot  are  a  resul?  of  the  diagrammatic  technique  and 
have  no  essential  significance.  The  particular  result  used  in  the  figure 
was  obtained  for  sample  HM45. 


It  is  most  noticeable  that  in  the  critical  regions,  where  equation 
(14)  is  obeyed,  the  branch  of  the  locus  for  T  >  T  is  separated  from  that 
for  T  <  T  ,  although  the  branches  are  parallel..  £his  only  arises  when  the 
crystal  structure  is  different  in  the  two  phases,  as  would  be  expected  for 
a  displacive  or  first  order  transition.  This  can  be  demonstrated  in' the 
following  maimer.  When  equation  (14)  is  obeyed  above  and  below  T  wTe  have 
that  c 


X<TC+) 


x*V) 


x(V> 


therefore 


-„«» 

-Up(T<)  ' 

U>  (T  **) 

-  p  c  ■ 


and  x(Tc“) 


Eyv’J' 


T  >  T 


(17a) 


(17b) 


which  become 


w’heo.  u  (T>) 


Hence  this  ratio  will  only  ha  unity  when  the  proportionality  constants 
in  equation  (16)  are  identical.  These  must  have  a  form  such  that  the 
susceptibilities  are  dimensionless  and  will  be  dependent  on  the  frequency, 
the  density  of  the  dipoles  and  their  dipole  moment,  among  other  features. 
Expression  (19)  has  been  checked  and  found  to  be  obeyed  for  these  samples. 
Figure  4  gives  the  ratio  of  relaxation  rates  in  the  two  phases  at  Tc« 

It  is  also  noticeable  that  the  small  grain  sized  sample  shows  a 
similiar  discontinuous  behaviour  to  the  other  samples,  although  close  to 
T  a  modification  appears  which  results  in  experimental  points  appearing  in 
tne  middle  of  the  two  parallel  branches.  A  similiar  effect  has  been  found 
for  individual  experimental  points  in  the  data  for  the  other  samples.  It 
can  be  generally  concluded  that  the  observed  dielectric  response  originates 
in  dopant  region^  undergoing  a  first  order  transition  which  is  forced  by 
the  host  lattice.  In  addition  to  this  the  dopant  regions  of  the  smallest 
grain  sized  material  will  assimilate  and  respond  together  with  the  host 
lattice  when  the  critical  slowing  down  of  the  host's  second  order  transition 
becomes  slower  than  that  of  the  dopant  regions  themselves. 

Outside  of  the  critical  region  relationship  (15)  has  been  shown  to 
apply  with  both  the  relaxation  rate  and  the  static  susceptibility  following 
activated  behaviour.  In  fact  it  can  be  seen  that  x(u-O)  takes  the  form 

x(u=0)  -  (T  +)  exp  {  -  nA/(kT)}  (20) 

c 

for  the  low  dopant  concentration  samples  (RM1A,  RM35A  and  RM35B) .  figure  5. 

The  frequency  index  n  arises  from  thg  gosling  of  the  relaxation  centre 
to  the  dipole  reversa^  modes  of  the  lattice  '  ’  and  measures  the  correlatio: 
between  these  motions  .  Outside  of  the  critical  region  the  centre  relaxes 
between  positions  separated  by  a  potential  barrier  creating  or  annihilating 
a  number  of  lattice  spin  reversals  in  the  process.  Such  local  spin  reversals 
determine  the  dipole  moment  of  the  centre,  but  notwithstanding  have  as 
independent  existence  as  propagating  excitations  on  a  time  scale  less  than  the 
relaxation  time  of  the  centres.  They  therefore  have  an  equilibrium  thermal 
population  which  determines  the  number  of  centres  that  are  excited,  and  which 
gives  equation  (15).  For  example  if  the  energy  per  dipole  reversal  is  kT 
then  A/kTc  spins  are  required  to  overcome  the  potential  barrier.  However  only 
a  fraction  n  of  these  are  retained  after  relaxation,  with  n  indexing  the 
coordination  of  the  centre  relaxation  with  the  lattice  motions. 


The  activation  energies  of  the  susceptibility  have  been  found  to  be 
2kTc  for  RM1A  and  R1I35A,  corresponding  to  a  single  ferroelectric  reversal  and 
3kTc  for  RM35B.  A  peak  of  5kT^  was  observed  for  the  0.67,  uranium  sample  and 
a  return  to  2kT£  for  the  maximum  concentration  of  1%,  figure  6.  The  values 
of  the  index  n,  as  a  function  of  concentration,  follow  a  similiar  pattern  vhicV 
is  given  in  figure  7.  Increases  in  centre  to  lattice  coupling  will  stiffen 
the  higher  frequency  lattice  vibrations  because  of  a  decrease  in  the  ability 
to  respond  on  their  time  scale.  It  is,  therefore,  to  be  expected  that  the 


ratio  of  relaxation  rates  below  and  above  T  will  also  follow  the  same  pattern 

c  r 

as  can  be  seen  in  figure  4. 


As  the  critical  temperature  is  approached  the  lattice  dipole  reversals 
couple  st^gngly  together  with  a  correlation  length  leading  to  the  critical 
phenomena  .  At  T  all  reversals  within  the  region  affected  by  the  centre 
will  take  part  and  the  susceptibility  x(T  +)  will  measure  their  number.  Since 
the  non-critical  value  of  equation  (15)  measures  the  thermal  population  as  a 
fraction  of  the  same  fixed  number  equation  (20)  results. 


Increasing  uranium  concentration  progressively  damages  the  lattice  in 
the  vicinity  of  the  uranium  centres,  as  can  be  seen  from  the  decrease  in  the 
activation  energy  of  the  relaxation  rate,  figure  8,  and  from  the  drop  in  T 
at  the  highest  uranium  concentrations.  This  is  also  observed  in  the  change  in 
the  relaxation  rate  ratios  shown  in  figure  4.  It  thus  appears  that  seme 
dopant  clustering  occurs  at  these  highest  concentrations.  At  0.6%  uranium 
(RM45)  this  clustering  is  insufficient  to  affect  the  host  lattice  ferro¬ 
electric  couplings  which  determine  T  and  the  uranium  centres  therefore 
involve  a  large  number  of  dipole  reversals  in  their  relaxation,  giving  a 
larger  value  of  n.  This  stronger  coupling  to  the  host  lattice  stiffens  the 
lattice  irself  and  increases  the  ratio  of  relaxation  frequencies.  At  1% 
dopant  (RM46)  the  uranium  centres  start  to  reduce  T  and  hence  must  be  damaging 
the  lattice  to  some  extent.  Since  u  (T  -)  is  much  5ess  than  to  (T  +)  it  appears 
that  this  damage  reduces  the  lattice^strf fness  through  the  proEucEion  of 
vacancies.  It  is  noticeable  here  that  equation  (20)  is  obeyed  if  ;<(T^+)  is 
replaced  by  x(T  ")  as  if  the  total  number  of  dipole  reversals  available  within 
the  same  sized  correlation  region  were  different  in  the  ferroelectric  and 
paraelectric  phases,  possibly  as  a  result  of  the  low  temperature  lattice 
distortion. 


Conclusions 

These  results  reveal  the  phase  transition  for  the  uranium  doped  ceramics 
to  be  generally  of  first  order.  It  is  possible  that  this  is  caused  by  eiectr- 
striction  from  the  undamped  grain  boundaries  modifying  the  second  order 
displacive  transition  expected  of  the  perovskite  structure.  The  sample  with 
the  smallest  grain  size  has  been  shown  to  revert  to  the  second  order.  The 
detail  of  the  main  ferroelectric  dispersion  however  could  not  be  observed  as 
it  lay  in  a  frequency  region  that  was  outwith  the  frequency  range  investigated. 


The  transition  characteristics  of  the  smallest  grain  sized  material 
indicated  the  presence  of  a  size  effect  in  limiting  the  approach  of  the  soft 
phonon  to  zero  frequency  at  the  critical  temperature.  An  estimate  of  the 
effective  grain  size,  and  hence  of  the  maximum  correlation  length,  can  be  made 
by  taking  th^  increase  in  the  permittivity  at  10  Hz  in  the  critical  temperature 
region  as  10°  and  using  the  observed  mean  field  values  of  one  half  for  v  and 
zero  for  n  in  equation  (2) .  A  value  of  thirty  to  forty  lattice  spacings  is 
obtained,  which  is  about  one  thirtieth  of  the  grain  radius  in  this  fine  grained 
material. 
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A  dielectric  response  originating  from  the  uranium  and  the  doj-ant 
sites  has  also  been  observed.  It  has  been  shovn  that  this  response  vas 
first  order  in  all  cases.  The  existence  of  dynamig  scaling  relationships 
was  verified  and  a  new,  non-critical,  relationship"  has  been  discussed  in 
qualitative  terms.  The  uranium  centres  have  been  found  to  couple  locally 
to  dipole  reversals  in  the  lattice.  Their  maximum  response  at  T  appears  co 
be  determined  by  the  number  of  dipole  reversals  available  for  coupling,  that 
is  by  their  domain  size  which,  except  for  the  0.4%  samples,  appears  to  be  of 
the  order  of  ten  lattice  spacings,  when  due  allowance  is  made  for  their 
concentration. 


The  single  calcined  material  represents  a  special  case  in  which  the 
uranium  domain  undergoes  a  first  order  transition  but  at  temperatures  close 
to  T  it  relaxes  faster  than  a  critical  slowing  down  second  order  process 
whiefi  it  then  follows  exactly.  In  this  way  it  becomes  part  of  the  host 
lattice  with  a  very  large  maximum^response  at  T  which  is  of  the  same  order 
as  that  of  the  permittivity  at  10  Hz.  It  is  therefore  very  likely  that  the 
critical  slowing  down  is  an  overdamped  soft  mode. 


In  summary  it  has  been  found  that 

r  * 

The  first  order/second  order  classification  is  a  consequence  of  a 
size  effect  and  due  to  the  volume /surface  ratio  of  the  grains.  This  has 
serious  consequences  for  the  mode  structure  of  the  material. 

There  is  no  direct  evidence  of  a  size  effect  in  the  correlated 
properties  of  the  materials  investigated.  In  all  cases  the  correlation  length 
was  a  fraction  of  the  observed  grain  sizes.  The  correlation  length  has  beer, 
estimated  as  being  approximately  thirty  lattice  spacings,  that  is  of  the  order 
of. one  hindred  angstroms. 
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Future  Work 


This  report  completes  the  work  that  will  be  carried  out  on 
ferroelectric  systems. 

During  the  remaining  six  month  period  of  the  present  contract 
it  is  intended  to  investigate  correlated  conduction  processes  in 
semiconductor  materials.  In  particular  silicon  will  be  considered  as 
the  model  semiconductor  and  the  interaction  of  correlation  size 
effects  with  carrier  relaxation  times  and  carrier  noise  will  be 
examined.  In  the  first  instance  the  investigation  will  be  theoretical  1 
based  but  experimental  data,  where  available  in  the  published  literatur 
will  be  made  use  of. 


Table  I 


Sample  Number  Mean  Grain  Diameter  Dopant  Concentration 


BM1A 

11  m 

nil 

EM35A 

2.5  m 

0.4% 

RM35B 

22  m 

0.4  % 

RM45 

24  m 

0.6% 

RM46 

- 

17  m 

1.0% 

Table  II 

Critical  exponents  for  the  ferroelectric,  high  frequency,  process 

below  and  above  T 

c 

e  -  (T  T)"y  e  -  (T  -  T  )"Y' 


Sample  Number 


Legends 


Figure  1.  Arrhenius  plot  of  the  zero  frequency  susceptibility  c£  the 
ferroelectric,  high  frequency,  process.  The  discontinuities  characteristic 
of  first  order  transitions  are  indicated  by  the  dotted  vertical  lines. 
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Figure  2.  Plots  of  the  permittivity  at  10  Hz  and  the  magnitude  of  the  low 
frequency  response  close  to  T  as  a  function  of  the  concentration  of  the 
uranium  dopant 

e  *T>T;oT<T  :  A(T)  +T>T;QT<T 

°»  c  c  c  c 

Figure  3.  Temperature  datum  trace  for  sample  R1145  showing  the  critical 
branches  of  gradient  -0.65  and  the  non-critical  branches  of  gradient  0.37. 
The  correlation  index  n  of  this  sample  was  measured  as  0.35  from  a  frequency 
plot  of  the  susceptibility.  The  direction  of  temperature  increase  is  shown 
by  the  arrow  heads. 


Figure  4.  Relaxation  rate  ratio  as  a  function  of  uranium  concentration 

Figure  5.  Arrhenius  plot  for  the  inverse  susceptibility  of  sample  BM35A. 
Note  the  infinite  temperature  extrapolation  of  the  activated  regions. 


Figure  6.  The  correlation  index  n  as  a  function  of  dopant  concentration. 

H  single  calcined  sample;  x  double  calcined  samples 


Figure  7.  The  activation  energy  of  the  susceptibility  as  a  function  of 
uranium  concentration 


Figure  8.  The  magnitude  of  the  activation  energy  of  the  relaxation  rate 
as  a  function  of  the  uranium  concentration. 
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